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Abstract 

In this paper, we investigate the speed of convergence and higher-order 
asymptotics of solutions to the porous medium equation posed in R iV . Ap¬ 
plying a nonlinear change of variables, we rewrite the equation as a diffusion 
on a fixed domain with quadratic nonlinearity. The degeneracy is cured by 
viewing the dynamics on a hypocycloidic manifold. It is in this framework 
that we can prove a differentiable dependency of solutions on the initial data, 
and thus, dynamical systems methods are applicable. Our main result is 
the construction of invariant manifolds in the phase space of solutions which 
are tangent at the origin to the eigenspaces of the linearized equation. We 
show how these invariant manifolds can be used to extract information on 
the higher-order long-time asymptotic expansions of solutions to the porous 
medium equation. 


1 Introduction 

1.1 Motivation 

The long-time asymptotic behavior of solutions is of fundamental importance in 
the study of nonlinear diffusion processes. Besides giving insight into characteristic 
qualitative properties, the long-time limit often sets a benchmark for the solutions’ 
expected regularity. The prototype for nonlinear diffusions is the porous medium 
equation 

d t u-Au m = 0 (1.1) 

posed in R A . Here m is a constant that we assume to be positive for a moment. As 
the name suggests, this equation in well-known for modelling the flow of gas through 
a porous medium, but also other applications are of relevance, e.g., as models for 
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groundwater infiltration, population dynamics or heat radiation in plasmas, cf. m 
Chapter2]. Solutions of (1.1) are assumed to be nonnegative. 

The qualitative behavior of solutions in general and the long-time behavior in par¬ 
ticular differs for different values of m. For instance, for m > 1, the diffusion flux 
mu m_1 vanishes where u = 0, and thus, if the initial configuration is compactly 
supported, the solution retains a compact support for all times. Hence, in such a 
situation, the porous medium equation features a free boundary value problem and 
the mass spreads slowly over the full space. We speak of “slow diffusion”. The 
situation is completely different if m < 1. Sticking to the mass-preserving range 
4^- < m < 1, solutions to (1.1) become positive instantaneously and remain posi¬ 
tive. In this range, (1.1) is called the “fast diffusion equation”. 

The characteristic long-time behavior of solutions to the porous medium equation 
(1.1) was first identified by Kamin in one dimension (29j [30] and by Friedman, 


Kamin, and Vazquez in several dimensions 1251 131] : Solutions approach the self¬ 
similar Barenblatt profiles [52, 7| 43] as they dwindle away to nothing. The optimal 
rate of convergence crucially depends on the choice of the initial datum. There are 
no uniform rates available without imposing a decay estimate on the initial datum 
at infinity, see |49, Theorem 1.3]. In the framework of entropy solutions, sharp rates 
were first found by Carrillo and Toscani [16], Otto [32] and del Pino and Dolbeault 

m. 


In the fast diffusion range, optimal rates of convergence and higher order asymp¬ 
totic corrections have been the subject of recent interest. In [2D], Denzler, Koch 
and McCann compute a number of higher-order modes of the long-time asymptotic 
expansions building up on a spectral calculation of Denzler and McCann [21] for 
the linearized dynamics. An overview on preliminary results is also provided in |20j . 
In the slow diffusion range, m > 1, the only known results on higher order asymp¬ 
totics are a full asymptotic expansion of Angenent [Tj in one space dimension, using 
the spectral calculation of Zel’dovich and Barenblatt [51], and an improved rate of 
convergence after centering the data by Vazquez [38], which is also proved in the 
one-dimensional setting only. Regarding the multidimensional dynamics, a prelimi¬ 
nary study was conducted by the author, who diagonalized the operator associated 
to linearized porous medium equation m- 

In the present paper, we successfully relate the information from [37] to the non¬ 
linear slow diffusion dynamics. This is achieved by exploiting ideas that were origi¬ 
nally developed for finite-dimensional dynamical systems. We prove that there exist 
finite-dimensional invariant manifolds in the phase space of the solutions, which are 
tangent at the origin to the eigenspaces of the linear operator, and every solution 
approaches these manifolds at a precise rate. In other words, we prove to what ex¬ 
tent the porous medium dynamics can be characterized by the dynamics of a system 
of ordinary differential equations. 

Invariant manifold theorems can be used to study the higher-order asymptotics for 
the porous medium equation: Once an order of decay is specified, our result ensures 
the existence of a finite-dimensional invariant manifold such that all solutions that 
are sufficiently close to the self-similar Barenblatt solutions approach that manifold 
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with at least this rate. Thanks to the discreteness of the spectrum, every eigenmode 
is accessible, and as a consequence, one can in principle compute the long-time 
asymptotics to any order. We will give four explicit applications in Theorems & 
PI below. 

The difficulty in making the dynamical systems approach work lies in the incompat¬ 
ibility between the space in which the spectral calculation was performed and the 
space in which solutions depend differentiably on their initial data. A framework 
in which the latter can be achieved was first developed by Koch in his habilitation 
thesis [36]. The main obstacle for establishing this differentiable dependency on the 
initial data is the variable support of solutions. In general, the domains on which 
solutions live differ from the one of the Barenblatt solution — the latter can be 
chosen stationary upon going over to self-similar variables. To overcome this, Koch 
suggests a nonlinear change of dependent and independent variables which results in 
a reparametrization of the graphs of the solutions. In the new coordinates, solutions 
live on the unit ball and Barenblatts become constants. The parabolicity of the 
resulting equation can be restored by trading the Euclidean metric for the liypocy- 
cloid metric, that is, the metric on the unit ball which is obtained when measuring 
distances along hypocycloids. It turns out that the solutions of the new equation 
depend even analytically on the initial data. To retrieve the information from 03. 
we have to modify the problem slightly for solutions which are far away from the 
constant attractor. 

The situation considered in the present paper differs from the one studied in [20] 
in many regards. Because solutions to the fast diffusion equation are positive ev¬ 
erywhere, it suffices to consider the evolution in relative variables to establish dif¬ 
ferentiable dependency on the initial data. Parabolicity is restored by carrying out 
the analysis on the asymptotically cylindrical cigar manifold, and the wcll-posedness 
theory for the rescaled equation is fairly standard. The main difficulty in [20] is the 
occurrence of continuous spectrum found in ED. above which only a finite number 
of eigenvalues are accessible. The authors need to work in weighted spaces to pene¬ 
trate into the spectrum. The choice of norms has a severe drawback: The resulting 
Banach spaces are not Hilbert spaces and the linearized operators are no longer 
self-adjoint. 

The present work can be considered as a model study for higher order asymptotic 
expansions for slow diffusions, such as the thin him equation with linear mobility 

d t u + V • (mVAm) = 0. 

The global attractor is again a Barenblatt solution HU m E3 HD]. A formal long 
time asymptotic expansion was obtained by Bernoff and Witelski [T2| in one space 
dimension, but no rigorous results on the higher order asymptotics are known. As a 
first step towards a rigorous understanding, building up on the spectral calculation 
in m, McCann and the author diagonalized the linear operator by exploiting the 
simple polynomial relation between the linear operators of the thin him equation and 
of the porous medium equation [31] . In a slightly simpler context, John constructed 
solutions to the thin him equation that depend differentiably on the initial data [2Sj • 
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1.2 Examples for higher-order asymptotics 


The applications of the invariant manifold approach, Theorems |l.lfjl.4| below, are 
most conveniently stated for the pressure variable in self-similar coordinates. By 
virtue of a standard rescaling argument, see, e.g., m Section 1.1], the porous 
medium equation can be equivalently written as a nonlinear Fokker-Planck equation 


d t u — A u m — V • (xu) = 0, 


( 1 . 2 ) 


which is often referred to as the “confined” porous medium equation. Under the 
same rescaling, the Barenblatt solutions w* (R, ■ ) become stationary, 

= M 2 ) + , 

m — 1 2 

where (/)+ = max{/, 0}, and the radius R of its support is determined by the mass 
of the solution. For studying the regularity of solutions, it is convenient to write 
the equation for the so-called pressure variable v = Upon changing the time 

scale from t to (m — l)f and setting a = the equation for v reads 

d t v — vAv — (cr + 1) (|Vu| 2 + x ■ Vu) — Nv = 0. (1.3) 

Existence and uniqueness of solutions to the Cauchy problem for the porous medium 
equation is established in [TO] - In [15.J, Caffarelli, Vazquez and Wolanski show that 
nonnegative solutions with compactly supported initial value have Lipschitz pres¬ 
sures after a certain time. In view of the Barenblatt pressure, this result is optimal. 
However, due to the existence of focusing solutions pm, one cannot hope to have 
Lipschitz regularity at any small time. 

To simplify the following discussion, we additionally normalize the “mass” of the 
initial configuration Vo such that 

J v ° +1 dx = j < +1 dx, (1.4) 


where n* = u*(l, ■). Then the long-time asymptotics of the porous medium equation 
is governed by v*. More precisely, 


lim || v(t) - v* 

ifoo 


= 0 . 


(1.5) 


We refer to Vazquez’ survey for a proof of (1.5). 


In this paper, we are interested in the characteristic behavior of solutions for large 
times, or equivalently, close to the Barenblatt solution. We thus select our initial 
data in a neighborhood of ry. More precisely, we shall assume that Vq is Lipschitz 
with 

N - p\\ L °°(v(v 0 )) < So and ||V(u 0 - p)\\l°°(t(v 0 )) < £o (1-6) 

for some <5 0 ; £ o > 0, and where, p(x) = — |x| 2 ). Regarding the validity of 

our main theorems, we could have equally chosen over p in (1.6). The actual 
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Notice that (1.6) entails that 


formulation, however, slightly simplifies the proofs. 

spt(uo) C 5 i +5 (0). 

Restricting ourself to the multidimensional setting (the one-dimensional case is con¬ 
sidered in HD, the spectrum of the linearized (and suitably transformed) equation 
consists purely of the eigenvalues 

\ ek = ( a + 1)(£ + 2k) + k{2k + 21 +N -2), 


where (£,k) G Nq x Nq, see Proposition 6.1 below. We relabel the eigenvalues in 


increasing order, that is, {A«,}h,fc)eN 0 xN () = {^j}jeN 0 with Aj < Aj+ 1- 

The knowledge of the spectrum (and the corresponding eigenvalues) of the linearized 


equation yields information on the higher-order asymptotics for (1.3) via the Invari 


ant Manifold Theorem 2.2 below: If j G Nq and A G (A j, A J+ i) are arbitrarily fixed, 


then the long-time asymptotic behavior of solutions up to terms of order e Xjt 
given by a finite system of ODEs, which result from restricting a transformed version 


is 


of (1.3) to a finite-dimensional invariant manifold. This manifold is tangent at the 
origin to the eigenspaces of the first j eigenvalues. In principle, we can compute 
the long-time asymptotics to any order. In the following, we make four explicit ap¬ 
plications: We derive the exact rate of convergence to the Barenblatt solution and 
compute the corrections in the long-time asymptotic expansions up to third order 
by restricting the dynamics to the invariant (sub-)manifolds. 

Our first result is a stability result for the Barenblatt pressure. 

Theorem 1.1 (Stability of the Barenblatt pressure). There exists > 0 with 


the following property. If Vq G O 0,1 satisfies (1.4) and (1.6) and v is the solution to 
(1.3) with initial value v 0 , then 


||u(f) — u*||l°°(r jv ) ^ e ( CT+1 b f or a u t > 0. 

Moreover, for any k G N 0 and [5 G N^ ; it holds 

II {v(t) - p ) || L ~(v(v(t))) < e~ (a+1)t for all t > 1. 

The same decay behavior in any spatial C k norm was already proved in [39j building 
up on Koch’s smoothness results [36] . 

The rate of converge in Theorem |l.l| is sharp. It is saturated by spatial translations of 
the rescaled Barenblatt pressure ty, that is, for any b G R jV , v(t, x ) = v*(x—e~^ +1 ^b) 


defines a solution to (1.3) which approaches u* with rate a + 1. 


The value Aio = cr + 1 is the smallest nonzero eigenvalue of the linearized operator, 
see Proposition |6.1| below. The corresponding eigenfunctions generate spatial trans¬ 
lations, which leave the porous medium equation (1.1) invariant. The eigenvalue 
Aoo = 0 corresponds to rescaling of mass, which leaves (1.1) invariant. Once the 


mass is fixed, see (1.4), this eigenvalues drops out of the spectrum 


The exact rate in Theorem 1.1 does not follow immediately from our invariant 


manifold theorem, Theorem 2.2 below, but requires a careful analysis of the nonlinear 


terms in the transformed equation, (2.5) below. 
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We obtain better rates of convergence when taking into account appropriate trans¬ 
lations of the Barenblatt solution. In this way, we get rid of the eigenvalue a + 1 
in the spectrum of the linearized operator, and the decay rate is given by the next 
eigenvalue 2(a + 1). 


Theorem 1.2 (First order corrections — spatial translations). Let A G (cr + l,2(<r + 
1)) be given. There exist £q, £o > 0 with the following property: If vq G C 0,1 satisfies 
( 1.4[ ) and (1.6) and v is the solution to (1.3) with initial value vq, then there exists 
a vector b G R A such that 


II v(t) - u*( • - e (<t+1) ^)||l-(r^) < e xt for all t > 0. 
Moreover, there exists a constant C > 0 such that 

Ri_ r( i)(0) C V(v(t)) C B 1+r{t) (0) for all t > 0, 
where r{t ) = e~^ <J+r>t \b\ + Ce~ xt . 


This result is obtained by a straightforward application of the invariant manifolds 
constructed in Theorem |2. 2 below, see Theorem 2.4 below, and a change of variables. 
For an analogous earlier result in one space dimension, we refer to |48j . 

The rate 2(<r+l), which it is not exactly achieved within this framework, corresponds 
to affine transformations of the Barenblatt solution. This is the first nontrivial 
eigenvalue as it does not reflect a symmetry of (1.1). Again, choosing appropriate 
affine transformations and centering the solution at the origin, we improve the rate 
of convergence once more. 


Theorem 1.3 (Second order corrections — affine transformations). Let A G (2(a + 
1), min{2(cr+l)+A^, 3(cr+l)}) be given. Suppose that a > — \. There exist £ 0 , <5 0 > 0 
with the following property: If v 0 G C 0,1 satisfies (1.4), (1.6) and 


/ 


Vq(xY +1 x dx 


0, 


(1.7) 


and v is the solution to (1.3) with initial value vq, then there exists a symmetric and 
trace-free matrix A G R v x R jV such that 


||u(t) — v*(x — Ax)\\l°° < e~ xt for all t > 0. 

Moreover, there exists a constant C > 0 such that 

B\— r (t') (0) C V{v(t)) C Ri+ r (p(0) for all t > 0, 
where r(t ) = e^ 2 ^ CT+1 ' )t |A| + Ce~ xt . 


Hence, the rate of convergence is given by the next eigenvalue in line, namely 
min{2(cr + 1) + N, 3(<r + 1)}. That is, we arrive at the first crossing of eigenvalues: 
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2(cr +1) + N corresponds to dilations or spatial rescaling of the Barenblatt solutions, 
and 3(a + 1) comes from pear-shaped transformations, see Figure 1 in pf7]. 

Our approach allows for a geometric interpretation of (1.4) and (1.7): By fixing 


the mass and centering the solution at the origin, we achieve that the first two 
eigenmodes are suppressed. In other words, the dynamics live on a submanifold of 
the invariant manifold which is parallel to the third eigenspace. 

The restriction to a > — \ is an artifact of the formulation of the invariant manifold 
theorem, and we will not make an attempt at getting rid of it in the present paper. 
In our final application of the Invariant Manifold Theorem |2.2[ affine transformations 
are modded out and we focus on dilations. 

Theorem 1.4 (Third order corrections — dilations). Suppose that N < a + 1 Let 
A G (2(<7 + l) + lV, 3(<7 + l)) be given. There exist e 0 , S 0 > with the following property: 


If v 0 G (F 0,1 satisfies (1.4), (1.6) and (1.7) and 


vo(x) a+1 x ■ Mxdx = 0 


( 1 . 8 ) 


for every trace-free matrix M G H NxN , and v is the solution to (1.3) with initial 
value vq, then there exists a constant c G R such that 

IK*) - R(t) 2_ “^ IU°° ~ e ~ xt f° r al11 ^ 

where a = 2 (a+i)+N an< ^ -^(*) = (l — ce _ ( 2 (°' +1 ) +Ar ^) a . Moreover, there exists a 
constant C > 0 such that 

Ri_ r(t )(0) C V{v(t)) C Ri+r(t)(0) for all t > 0, 

where r(t ) = e _< ' 2( - <T+1 ^ +iV - ) *Q:|c| + Ce~ xt . 

In principle, similar statements on corrections of even higher order could be derived 


via the Invariant Manifold Theorem 2.2, but the computations would get more and 
more involved the larger the eigenvalues were chosen. 

The remainder of the paper is organized as follows: In the following section, we 


express the rescaled pressure equation (1.3) in new variables that are suitable for 


our dynamical systems approach to work. We present a theorem on wcll-posedness 
of the resulting equation and state our main theorem about invariant manifolds. 
Finally, we reformulate Theorems |1.1]-|1.4| in terms of the new variables. Section 


[3] contains the proofs of Theorems 1.1-1.4 as a consequence of the corresponding 
statements for the new variables. The rest of the paper is entirely devoted to the 
analysis of the transformed equation. In Section [4| we study its linear version and 
establish a well-posedness theory that is appropriate for addressing the nonlinear 
equation with a fixed point argument. This is the content of Section [5] The invariant 
manifold theorem is proved in Section [6j and the applications are studied in Section 
□ 
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2 New variables and main theorem 


The main difficulty of applying dynamical systems arguments to the porous medium 


equation (1.1) (or (1.3)) consists of finding a setting in which solutions depend dif- 
ferentiably on the initial data. One problem appears to be the free moving boundary 
at which the regularity of solutions breaks down and at which perturbations around 
the (stationary) Barenblatt solutions are hard to define. The elegant mass-transport 
approach chosen in [21] 07j and which builds up on Otto’s gradient flow interpreta¬ 
tion [42] overcomes this obstacle via Brenier’s theorem mm, but a differentiable 
dependency on the initial data using mass transport techniques seems out of reach. 
A framework in which the latter can be achieved was introduced by Koch in his 
habilitation thesis [36]. Koch suggests a nonlinear change of dependent and inde¬ 


pendent variables which transforms the free boundary problem (1.3) into a partial 


differential equation on a fixed domain. Moreover, in the new variables, Barenblatt 
solutions (A, •) are mapped onto the constant R. 


With total mass fixed as in (1.4) (i.e., R — 1), at any time t > 0, the transformation 
of the spatial coordinates reads 


x 


z = 


a/2 v(t, x ) + |x| 2 


( 2 . 1 ) 


It is clear that the new coordinate z lies in the unit ball, and the transformation 
reduces to z — x if v is the Barenblatt solution u*. Notice that ( z , y / 2n* ( 2 )) is 
obtained by projecting (x, \j2v[t, x)) orthogonally onto the graph of y/2v*, which is 
a half-sphere. 

We introduce the new dependent variables as perturbations around the Barenblatt 
solution. More precisely, if w(t,z ) denotes the distance between (x, \J2v(t, x)) and 
the half-sphere, it holds: 


1 + w(t, z) = a/2 v(t, x) + |x| 2 . 


( 2 . 2 ) 


As announced above, this change of variables transforms the Barenblatt solutions u* 
into the constant function 1, or, in other words, perturbations are constantly zero, 


i.e., w = 0. The change of variables will be made rigorous in Lemmas 3T and 332 
below. These two lemmas in particular show that there is a diffeomorphism between 
the graph of v and the graph of w. I 11 Figure [l] above, we illustrate this change of 
variables around the Barenblatt solution. 

For further references, we note that the inverse change of variables reads 


x 


= (w(t,z) + l)z, 


v(t,x)—p(x) = w(t,z) + -w(t,z) 2 . 
The transformed equation for the to-variables is 
d t w - p~ a V • (p a+1 Vw) 


(2.3) 

(2.4) 


= (N + 2a + l)p 


|Vw/ 


1 + w + z ■ Vw 


- p~° V • 


p^z 


| Vw 2 


1 + w + z ■ Vw 


(2.5) 


8 

















Figure 1: The plot on the left shows the graph of the Barenblatt pressure (blue 
curve) and two perturbations in the (x, v) -plane. The plot on the right shows the 
corresponding graphs in the (z, 1 + ?n)-plane. 


see also Lemma 3.2 below. Since w represents the perturbation around the constant 


(former Barenblatt) solution, we will often call (2.5) the (nonlinear) perturbation 
equation. 

In [36], Koch develops a regularity theory for equations of this type in Holder spaces. 


For a slightly simpler version of (2.5) in which the nonlinearity is independent of 
the function itself (and thus only depends on derivatives), Kienzler [3T] improves 
on Koch’s results and establishes a theory for small Lipschitz initial data. Building 


up on Kienzler’s techniques, we will construct smooth solutions to (2.5) for small 
initial data in C 0,1 . These solutions are uniformly controlled in any C k norm by 
the Lipschitz norm of the initial datum, which prevents the nonlinear terms from 
degenerating. In addition, solutions depend analytically on the initial datum. More 
precisely, we prove 

Theorem 2.1 (Well-posedness for the perturbation equation). There exists e > 0 
and 5 > 0 such that for every g G B £t s, there exists a unique solution w to ( |2.5[ ) . This 
solution is smooth and depends analytically on g. Moreover, ||ta|| 

L°° + IMkip < i 

and 

i‘ + W|5*afv»(t,n|< WIlip. (2.6) 

for any (t , z ) G (0, oo) x -£>i(0) and k G N 0 and j3 G N, 


N 


Here and in the following, we use the notation 

B £ j ■■= { w G C 0,1 : \\w\\ L oo < S, 11 rail Lip < ej ■ 

The analyticity of the solution is a byproduct of the fixed point argument applied 
in the construction of solutions of the nonlinear equation. The underlying idea has 
been first used by Angenent m and has been further exploited by Koch and Lamm 


Before stating our results on invariant manifolds and long-time asymptotics, we 
need some preparations. Let £(£) = {Afc}fc g N 0 be the (discrete) spectrum of the 
linear operator C = —pA + (er + 1 )z ■ V on the Hilbert space L 2 a = L 2 (dp a ), where 
dp a = p a dx. See Proposition |6.1| below for more details on the spectrum of C. 


Suppose that the eigenvalues are in increasing order and not repeated, so that < 
Afc+i. For K G N 0 arbitrary, consider the subspace E c of L 2 that is spanned by 
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the eigenfunctions corresponding to the first K eigenvalues, and denote by E s its 
orthogonal complement, i.e., L 2 = E c ® E s . Let A and £ gap be positive constants 


so that A e (e < A +i) _|_ £ 


gap, 
t 


e k — e 


gap 


). The semi-flow associated to (2.5) will be 


denoted by S t , that is S t (g ) = w(t) if w denotes the solution of (2.5) with initial 
value g given by Theorem 2T above. 

We are now in the position to state our main result. 


Theorem 2.2 (Invariant manifolds). There exist constants e, 5 > 0 and e, 5 > 0, 
a function 6 : E c —>■ E s and manifolds W c = {g c + 9(g c ) : g c £ E c } and W^ oc := 
W c fi B e> s with the following properties: 

1. The function 9 is a differentiable map from B s .s to C 0,1 satisfying 9(0) = 0, 
D9( 0) = 0 and \\9(g c )\\ c ov < Us'dlco.i f or al1 g c £ E c D B £jS . 

2. If g £ W c fl B £ ~ s , then S t (g) £ W l £ c for all t > 0. Moreover, there exists a 
unique g £ W £ oc such that S t (g ) £ W l ° c for all t > 0 and 

11^)A* for all t> 0. 


3. If {w(t )} t <o is a negative semi-orbit of (2.5) with g £ W c and w(t) £ B £: s for 


all t <0, then w(t) £ W £ for all t < 0. 


The first property shows that the local center manifold W £ is tangent at the origin 

/■loc 

J C • 111C OCL/UllU piupci IJ D11WVVO tllCXL q 


to the eigenspace E c . The second property shows that W £ oc is locally flow invariant 


under the semi-flow {5'(t)} t >o for the nonlinear perturbation equation (5.1) and 


shows to what degree solutions to the perturbation equation (2.5) are characterized 
by a system of ODEs. 

We finally give four applications of the invariant manifold theorem. In the first 
example, we study the long-time limit of small data solutions of the perturbation 
equation. 

Theorem 2.3 (Stability). There exists positive constants e and 5 with the following 
properties: If g £ B £ ~ & and w solves (2.5) with initial datum g, then there exists a 


constant a £ R such that for all k £ Nq and (3 £ it holds 


\\dt d z (M*) - «) 




for all t > 1. 


(2.7) 


The rate of convergence in the above theorem is sharp. It is saturated by spatial 
translations. The exact rate does not follow automatically from Theorem 2.2, but 
requires a detailed analysis of the nonlinear term. 

Notice that any constant (except 1) is a solution to the perturbation equation (2.5). 
As a consequence, the constant long-time limit is not necessarily zero. On the level 
of the porous medium equation (1.1), this corresponds to the invariance of solutions 
under the rescaling of mass. Once the mass is fixed as in (1.4), this constant drops 
out in (2.7). In this way, we obtain the first order corrections in the long-time limit 

of (pi). 
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Theorem 2.4 (First order corrections — spatial translations). Let A G (cr + 1, 2(cr + 
1)) be given. There exists positive constants e and 5 with the following properties: 


If g G B~~ 5 and w solves (2.5) with initial datum g and satisfies 


lim / w(t) dp a = 0 , 
ttoo J 

then there exists a vector b G R A such that 

||tc(t) — e~^ c+v>t b • ^Hcfo.i < e~ xt for all t > 0 . 

The rate of convergence can be further improved when the center of mass is asymp¬ 
totically fixed at the origin. In this case, spatial translations are suppressed and 
the dynamics are governed by affine transformations. Thus, we climb the eigenvalue 
ladder one rung up: 

Theorem 2.5 (Second order corrections — affine transformations). Let A G (2(cr + 
1), min{ 2 (cr + 1 ) + N), 3 (cr + 1 )}) be given. There exists positive constants e and 5 


with the following properties: If g G B~~ 5 andw solves (2.5) with initial datum g and 
satisfies 


lim / w(t) dp a = 0 and 

tfoo J 


lim / w(t)zdfi a — 0 . 

ifoo J 


Then there exists a symmetric and trace-free matrix A G R NxAr such that 


IN*) 


_ -2(«r+l)t 


^ • Az\\ C o, 


< p~ M 


for all t > 1 . 


We have a fist eigenvalue crossing when iV = er + 1: It holds that A 3 = min{2(cr + 
1) + N,3(cr + 1)}. On the level of the porous medium equation, the eigenvalue 
2(a + 1) + N corresponds to dilations, which is, by the invariance of the porous 
medium equation under time shifts, an artifact of the dynamics. We can access 
the corresponding eigenvalue by moding out certain second order moments that 


are preserved under (1.1). On the level of the perturbation equation (2.5), this is 


achieved in the following 

Theorem 2.6 (Third order corrections — dilations). Suppose that N < a + 1. Let 
A G (2(cr + 1) + N, 3(cr + 1)) be given. There exists positive constants e and 5 with 


the following properties: If g G B~$ and w solves (2.5) with initial datum g and 
satisfies 

lim / w(t) dp, a- = 0 , 

ifoo J 

lime^ +1 ^ / w(t)zdp a — 0, 

t'j'OO J 

lime 2(<T+1)< J w(t)M:z® zdp a = 0 

for every trace-free matrix M G H NxN . Then there exists c G R such that 


| w{t) - e - 2 ( 0 +b +Ar )* c (i - 7 | 


I c° 


1 <e xt for all t > 1 , 


where 7 = N 1 (2(cr + 1) + N). 
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In principle, every eigenvalue is accessible through the Invariant Manifold Theorem 
m However, as the complexity of the corresponding eigenfunctions increases with 
every eigenvalue, we refrain from studying further examples. 
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3 Higher order asymptotics for the pressure vari¬ 
able 

In the present section, we show how the statements on the long-time asymptotics 


for the perturbation equation (2.5), Theorems 2.3 2.6 imply the corresponding 


statements on the long-time asymptotics for the confined pressure equation (1.3), 
Theorems 11.1111.d[ 

It is convenient to recall and introduce some notation. For any function v : [0, oo) x 
R N we denote by V{v) its positivity set and by T(y(t )) its positivity set at 

time t, that is, Viv) = {(i, x) G [0, oo) x R A : v(t,x) > 0} and V(y(t)) = {;r G 
R iV : v(t, x) > 0}. 


In a first step, we show that the transformation (2.1), (2.2) is well-defined and show 


how smallness conditions for Vo in (1.6) turn into smallness conditions for the new 
variables. 


Lemma 3.1 (Change of coordinates I). Suppose that vq G C 0,1 satisfies (1.6) for 
some £q, 5q > 0 with Eo + 2<5 0 < 1. Let <3>o : V{vq) —> Ri(0) be given by 


‘hnbC = 


x 


\/2v 0 (x) + |x| 2 ’ 
for all x G V(vq). Then the following holds: 

1. The mapping $o is a dijfeomorphism. 

2. Let g : T>i(0) —* R be given by 

g{®o(x)) = y/2v 0 (x) + |x| 2 - 1, 
for all x G V(v 0 ). Then 

1 + 25q 


M L°° < 25 0 and 


Lip < 


1 — 2 J 0 — 


£ 0 - 


Proof. 1. The Jacobi matrix of <l>o is given by 

1 x <g) (X7v 0 (x) + x) 


V$ n m = 


sj2v 0 (x) + |x| 2 (2v 0 (x) + |a;| 2 ) 3 ^ 2 


where 1 denotes the identity matrix in R A x R A . We compute its Jacobian deter¬ 
minant: 


= (2v 0 (x) + \x \ 2 ) 
= (2v 0 (x) + | a; | 2 ) 


x ® (Vu 0 (^) + x )\ 
2v 0 (x) + \x\ 2 ) 

PL_ 1 

2 (2v 0 {x) - x ■ Vv 0 (x )), 


N 

2 det 
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det V$o(^) 





















where the second identity follows from the auxiliary formula^ det (1 + a <E> b) = 
1+a-b. Because 2u 0 +| a; | J = 2(n 0 —p) + l and 2vq+x-Vvq = 2(vo—p)-x-V(vo—p) + l, 


we deduce from (1.6) that the Jacobian is positive on the positivity set of vq. Hence 
$o is a diffeomorphism. 

2. Notice that g is well-defined because <ho is a diffeomorphism. The uniform bound 
on g is a consequence of the first assumption on Vo in (1.6). For the Lipschitz bound, 
notice that 

1 + 9(z) 


Vv 0 (x) + x = 


1 + g(z)+z-Vg{z) 
if z — < ho(®)- Therefore, multiplying by x = (1 + g)z gives 


Vs(*) 


1 + 9 


1+g + z-Vg 


= 1 


x 


■ V Vn + \x\ 2 2v t 


X 


2v 0 + \x\ 2 


_Vno 

2v 0 + \x\ 2 


and the statement follows via (1.6) 


The following lemma justifies the inverse transformation (2.3), (2.4). Moreover, we 
show that under this transformation, solutions to the perturbation equation (2.5) 
become solutions to the confined pressure equation (1.3). Thanks to the uniqueness 
of the associated initial value problems, both equations are thus equivalent. 


Lemma 3.2 (Change of variables 11). Let e, S be as in Theorem 2.1. Suppose vj 
is the solution to (2.5) with initial value g G B £ j. Consider T : (0, oo) x T>i(0) —* 
(0, oo) x R a given by 

T(t, z) = (t, (1 + w(t, z))z), for (t, z) G (0, oo) x B i(0). 

Then the following holds: 

1. T is a diffeomorphism onto its image. 

2. Let v be defined by 

p(z)( 1 + w(t, z)) 2 if (t, x) = T(f, z) for some z G -Bi(O), 


v(t, x) = 


0 


otherwise. 


Then v is C 0,1 and solves the confined pressure equation (1.3). 


Proof. 1. We start with the computation of the Jacobian matrix of T, 


1 0 

d t (t,z)z z <8) Vw(f, z) + (1 + w(t, z))l 


1 This formula immediately follows from the calculation 


1 0 

b T 1 


1 + a ® b a 
0 1 


1 0 

1 


1 a 

0 1+a-b 
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Arguing similarly as in the proof of the previous lemma, we deduce that 


det(d t \I/(£,z), V\h(f, z)) = (1 + w(t,z )) (1 + w(t, z) + z ■ Vw(t, z )), 


and the expression on the right is positive by Theorem 2.1 Thus T is a diffeomor- 
phism onto is image. 

2. It is clear that a is well-defined, bounded and Lipschitz on R A because la is 


bounded and Lipschitz. To verify that v indeed solves (1.3), we give a quite formal 
argument. Starting point is the identity 


a/2 v(t,y) + |/ 2 = 1 + w It, 


x 


a/2 v(t,x) + \x\ 2 ) 

A multiple use of the elementary differentiation rules yields 

(1 + w + z- Via) d t v = (1 + w) 2 d t w , 
(1 + w + z ■ Via) (Va + x) = (l + ia)Via, 


and 


(1 + w + z ■ Vw) (Aa + N ) 

2z ■ ( D 2 w'Vw ) + | Via | 2 | Via | 2 (2z ■ Vw + z <S> z : D 2 w ) 


= A w — 
Noting that 

2 -V 


| Vrap 


1 + w + z ■ Via 


2z ■ ( D 2 w\7w ) 


(1 + w + 2 • Via/ 


| Via p 


1 + w + z-VwJ 1 + w + z-Vw 1 + w + z-Vw 
the previous formula simplifies to 

| Via | 2 


(2z ■ Via -\- z® z \ D 2 w ) , 


(1 + in + z • Via) (Aa + N) = Aw+(N— 1) 


1 + la + z ■ Via 


V- z 


| Viap 


1 + w + z ■ Via 


Multiplying both sides of the equation by (1 + la) 2 a = p(z), a short computation 
leads us to 

(1 + in )' 2 (1 + w + z- Via) (aAu + Nv) 


= V ■ ( p ( Via — ^ 


| Viap 


+ z ■ Via — \z I 


1 + la + z ■ Vw 
| V la | 2 


1 + w + z • Vw 


+ (N- 1 )p 


\Vwf 


1 + w + z ■ Via 


Without much effort we can now derive (1.3). 


By the virtue of the preceding two lemmas, we are now in the position to deduce 


Theorems 1.1 1.4 directly from Theorems 2.3 2.6 
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Proof of Theorem H3 Let e and 5 be as in Theorem 


2.3 


Let g be defined as in 


3.1 


Then g G is £q and <5 q are chosen sufficiently small. Let w be the 


Lemma 

solution to (|2.5|) with initial value g. By Theorem 2. 3[ there exists a constant a G R 
such that 


|| w(t) — < e for all t > 0,, 

\\d^(w(t) - a)\\ L « < for all f > 1,, 


(3.1) 

(3.2) 


for all k G N 0 and /3 G N/. In particular, w{t ) —> a for large times. Since on the 
other hand 

v — p = w + \w 2 inV(v), (3.3) 


by Lemma 3.2, and v —> c* = (p) + by (1.5), we must have that a 2 + a = 0, and thus 
a —— 0 since 11^11//°° < 6 < 1. It immediately follows from (|3. 11) and (|3.3|) that 


| v(t,x) — p{x) | < e ( CT+1)f for all x G V(v(t)), t> 0. 
This estimate in particular implies that 


(3.4) 


1 - < Id < 1 + Ce~ 


for all x € dV(v(t)) and t > 0 and some C 
V(v(t)) c B 1+ Ce-(<r+i)t (0) 


> 0. Therefore B 1 _ Ce -^+ i)t(0) C 


Moreover, in (3.4), we can easily replace p by ti* = 
(p) + because v — c* = v < v — p outside of -£>i(0). On the other hand, since 
-Bi(O) \ B(v(t)) C B\ (0) \ it holds that |c*(x)| < ou tside of 

V{v(t)). In any case 


| v(t,x) — n*(x)| < e O+L* for all t > 0. 

To deduce the corresponding statement for higher order derivatives, we argue as 
follows. From the proof of Lemma 322 we recall that 

1 + w 


di(v - p) = 


1 + w + z ■ Vic 


diW. 


The expression on the right-hand side and all of its derivatives with respect to z 
are bounded by e~( a+1 ') t thanks to (3.2). This implies the desired statement for 
(k, |/3|) = (0,1). Moreover, derivatives of the diffeomorphism 

. , s x 

z = $ t (x) = 


a/2 v(t, x) + \x\ 2 

are uniformly bounded, for instance, 

1 x <g> (Vc + x) 1 — z <E> Vic 


VMh = 


\/2c + 


x\ 


a/2i7+~ 


x\ 


1 + w 


16 


























and the control of higher order derivatives follows via the chain rule and (3.2). 
Hence, for any (3 G Nq\ it holds that 


v ~ p)h^(v(v(t)) < e {a+1)t for all t > 0 , 


and the control of the temporal derivatives follows upon using the equation (1.3) in 
the form 

d t v = vA(v — p) + (a + l)Vn • V(v — p). 


This concludes the proof of Theorem 1.1| 


Proof of Theorem \l/z\ We let § and 3 be as in Theorem 
Lemma 


2.4 


and define g as in 


3.1 


We suppose that £o and 5 o are small enough so that g G B s $. From the 

(3.5) 


statement of Theorem 241 and the proof of Theorem 1.1 , we deduce that 

I loo < e ~^ +V)t for all t > 0 , 


\w 


and thus, 


lim / w(t) dp a = 0 . 

t 'foo / 


Consequently, thanks to Theorem |2.4| there exists a vector h G R v such that 


\w 


(t) - e~ {rT+1)t b ■ z\\co,i < e~ xt for all t > 0 . 


(3.6) 


Thus, Lemma 3.2 and ( |3.5 ) imply that 

v{t,x) - p(x - e~ {a+1)t b) = w(t,z ) - e~^ +1)t b ■ z + 0(e~ xt ), 


for all (t,z) G V(v), and the right-hand side is of order e xt by (|3.6|). Finally, as in 


the proof of Theorem 0 we can easily switch from p to v* = (p)+ on the left-hand 
side, extent the statement to all x G R A and derive estimates on the positivity set 
of v. m 


Proof of Theorem 1.3. We start recalling the well-known fact that the confined 


porous medium equation ( 1 . 2 ) preserves zero-center of mass. Indeed, a short com¬ 


putation (performed on the level of (1.3)) reveals that 


-f- [ v a+1 Xi dx = 

dt ./ 


-0 + 1 ) 


,<T+1 


Xi dx, 


and thus, the centering condition (1.7) is valid for all times. Let e and 5 be as in 


Theorem 2.5, and define g as in Lemma |3.1[ We choose £q and <5 0 sm all enough so 


that g G B~~ 5 . Let A G (max {a + 1, l}2(cr + 1)). By Theorem 
vector b G R jV such that 


2.4 


there exists a 


ft,x) — vjx — e (<T+1)t 6) | < e xt for all x G R jV ,f > 0. 


(3.7) 
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We claim that 6 = 0. Indeed, by the symmetry of w*, it holds that 




*{x) a+1 (x + e (cr+1)t 6) dx 


\b\ ~ J v*(x) a+1 e~ i ' C7+1 ' lt bdx = 

By a change of variables and because v is centered at zero, we further have that 
e -o+i)/|kl < f \ v ^ x — e ~t (T+i ^ t by +l — v(t,xy +l \\x\dx. 


In the case a > 0, an application of (3.7) yields \b\ < e ( <T+1 ^ t , and thus, letting 


i t oo, we see that 6 = 0. Likewise, if — ~ < a < 1, we apply (3.7) and obtain 


| 6 | < e 1 K< j + 1 )* ; a nd again, 6 = 0 follows upon passing to the limit in t. Because 
6 = 0 , it holds that 

lime^ +1 ^ / w(t,z)zdfi a = 0, 
ttoc J 


and thus, Theorem 2.5 is applicable. The argument of how to deduce Theorem 1.3 


from Theorem 273 proceeds similarly as in the proof of the previous theorem. We 
omit the details. ■ 


Proof of Theorem l.f The proof is conceptionally very similar to the proof of the 


previous theorem. We notice first that the second moment condition in (1.8) is 


preserved under the evolution (1.3) because 


^ f v a+1 x ■ Mxdx = — 2(a + 1) f v a+1 x-Mxdx 


for every trace-free matrix Mg R 


NxN 


Let i, 6 be as in Theorem 


2.6 


and define 


3.1 


We choose e 0 and 5 0 small enough so that g G B i $. Let 


g as in Lemma 
A G (2(a + 1), 2(a + 1) + N). Hence Theorem 1.3 implies that 

| v(t, x) — v*(x — e~ 2( ' a+1 ' )t Ax)\ < e~' xt for all x G R iV , t > 0, 


(3,8) 


for some symmetric and trace-free matrix A G H NxN . We claim that A — 0. Indeed, 
we can find a trace-free matrix M with \M\ ~ 1 and 
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u*(x) CT+1 a; • ( MA)xdx 


Since u* satisfies (|1.8|) by symmetry, we thus have 

e -2(a+l)t\ A \ _ 


+ e~^ a+l)t Ax) ■ M(x + e ~ 2(a+l)t Ax) dx 


+ 0{e~^ +1)t ) 


<{x — e 2( ' c+v>t Ax) a+1 x ■ Mxdx 


+ 0{e~ 4{rT+1)t ). 


Invoking (1.8) for v, we further have 


e - 2(a+1)t \A\ < ||x;(^) CT+1 - v*(x - e~ 2 ^ +1)t Ax) a+1 \\ L oo + 0{e~ i{a+1)t ). 
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Then, because a > 0 by assumption, an application of (3.8) yields that e~ 2 ^ Jrl '^ t \A\ < 
e~ xt f or a jj ^ > o, and thus A > 2(cr + 1) implies A = 0. The remainder of the proof 
is very similar to the proofs of the previous theorems. We only remark that 


v(t, x ) - -(t - |x| 2 ) - e-^+d+AO^f _ j\x\ 2 ) = 0(e~ xt ) 


follows from Theorem 2.6 We then compute 


2-R(f) a 


r M«) : 


\X\ 


2 ) = ^(1 - \x\ 2 ) + C e-( 2 ( ff+1 ) +A T( 1 - 7 |a:| 2 ) + 0(e ~ xt ), 


where we have used that (1 — 2a) 7 = 1. The statement of the theorem now can be 
deduced as before. ■ 
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4 The linear problem 


In this section, we study the initial value problem for the degenerate parabolic 
equation 

d t w - p~ a V • (p CT+1 Vw) = / (4.1) 


in -£>i(0), where p(z) = |(1 — \z\ 2 ) and a = 


2 —m 


-l,oo). A big effort is made 


2 V 11/ m— 1 

to obtain suitable regularity estimates that serve as a basis for the well-posedness 
theory developed later in Section [4] for the nonlinear perturbation equation (2.5). 
Here, the choice of the topology is crucial. To reach the ultimate objective, the 
construction of invariant manifolds for the nonlinear equation in Section [6j a differ¬ 
entiable dependence of solutions on the initial data is necessary. For nonlinearities 
of the form 

f = pF- p~°V ■ (p" +1 zF) 

in which F is a quadratic function of Vw, we expect (and, in fact, will prove) 
that such a differentiable dependence can be established for solutions with small 
Lipschitz initial data. In order to obtain suitable regularization estimates, we need 
to work with space-time norms that deteriorate as t f 0. We thus naturally arrive 
at Carleson measures. 

The analysis of the present section is inspired by Koch and Lamm’s recent ap¬ 
proach to tackle semilinear parabolic equations with rough coefficients [37] and, 
even more, by Kienzler’s adaption of this approach to a subelliptic parabolic equa¬ 
tion [33], which closely resembles ours. The authors use standard tools from the 
theory of constant coefficient linear equations such as (Gaussian) decay estimates 
and Calderon-Zygmund-type estimates to establish bounds on the solutions in cer¬ 
tain Carleson measures that reflect the regularity theory for the corresponding linear 
equation. The underlying idea goes back to Koch and Tataru’s work [38], where a 
Carleson measure formulation of a BMO norm turned out to be the crucial ingre¬ 
dient in the study of regularity for the small datum Navier-Stokes equation. Many 
ingredients from [33] appeared earlier in Koch’s habilitation thesis [3B], where Koch 
established a Schauder theory for Holder continuous initial data. 


The main difference between (4.1) (and likewise the equations studied in [331 53] ) 


and the equations considered in [37] is the failure of strict parabolicity of (4.1) at 
the boundary of the domain. However, we can easily overcome this problem by 
the following observation: Equation (4.1) can be interpreted as a heat flow on a 


certain so-called weighted manifold, that is, a Riemannian manifold to which a new 
volume element is assigned, typically a positive multiple of the one induced by the 
Riemannian metric. The theories of weighted manifolds and heat flows thereon can 
be developed parallel to the Riemannian counterparts, see, e.g., [SB]. In particular, 
Gaussian decay estimates on the heat kernel exist and a Calderon-Zygmund theory 
is available. We do thus expect that suitable Carleson measure estimates can be 
established if we carry over Koch and Lamm’s approach from the Euclidean to the 
Riemannian setting. In fact, we will still work on the Euclidean unit ball, but we will 
alter the metric accordingly. By working with the new metric, which is a Carnot- 
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Caratheodory distance, see, e.g., [9], we restore the parabolicity of the equation. 
In the context of the porous medium equation, the application of such Carnot- 
Caratheodory distances has been proved useful in the pioneering works [18], [36], and 
recently in [3lj. See also E23 and (28] for a similar perspective on the fast diffusion 
equation and thin film equation, respectively. 

Our goal is the following result: 


Theorem 4.1. Let p > max {iV + 2, For every function g G C 0,1 and f G 

Y(p), there exists a unique solution w to (4.1). Moreover, w G X(p) and 


' w l|v:(p) + II^II Lip ll/ll>'(p) T \\g\\ Lip* 


Here, we have used the notation X(p) = {u> : ||tc||x(p) < oo} and Y (p) — {/ : 
\\f\\Y(p) < oo}, where 

IMUfr) := sup \Q d r {z)\~p (r 2 \\dtXw\\ LPiQ d Az)) + r(r + \Z~p[z))\\X 2 w\\ LP{Q d r{z)) 

^eBi(o) v 

0 <r<s/2 

+ ^ 2 ||pV 3 u;|| i p( (9 d( 2 ))j 

+ sup (||^V-u;|| L p(g( T )) + ||V 2 u;|| L p(q( T )) + \\pX 3 w\\ LP (Q( T ))) , 

T> 1 

\\f\W(p) ■= sup \Q d r (z)\-v (-^==||/|| L p (Q d ( 2 )) +r 2 ||V/|| L p (Q d (2)) ) 

^eSi(o) \ r + V P\ z ) ) 

0 <r<V^ 

+ sup (||/||lp(Q(T)) + I|V/||lp(Q(T))) , 

T> 1 


and where Q^(z) and Q{T ) denotes the (intrinsic) time-space cylinders, 


Q d r (z ) = ,r 2 J x B d r {z) and Q(T ) = (T,T + 1 ) x ^(0). 

Here, the intrinsic balls are relatively open subsets of the closed Euclidean unit 
ball, B d (z) = < z' E B\{z ) : d(z,z') < r > and d is the new metric on I?i(0) which 


is induced by the heat flow interpretation. In particular, B d (z) = £>i(0) for every 
sufficiently large r. Finally, the Lipschitz norms are taken with respect to the spatial 
Euclidean topology, that is, ||w||Li P = || Vw\\l°c,. 

By the linearity of the equation, we may and do split the problem into the ho¬ 


mogeneous problem, where / = 0, Proposition 4.2 and the problem with g = 0, 


Proposition 4.4 From these two cases, Theorem 4.1 follows by superposition. 


Though an abstract theory for weighted manifolds is well developed, cf. |2S] and 
references therein, we will present a self-contained theory in the subsequent sub¬ 


sections since, with regard to the nonlinear equation (2.5) to which the regularity 


estimates of this section will be applied, results have to be tailored to our needs. 
We conclude this introduction with a precise definition of weighted manifolds and the 


description of the heat flow interpretation of (4.1). Suppose (PV4, g) is a Riemannian 
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manifold and uj a positive function on A4. If we denote by dvol g the volume element 
on (A4,g) and set dji = cudvolg, then the triple (A4,g,/u) constitutes a weighted 
Riemannian manifold. If the Laplacian is defined via the Dirichlet form on (A4 , g, p), 
its representation in local coordinates reads 


N 


A^w = 


"V 3if s£i 


V ft (g y u\/detgft 


w 


where {g* J }p- is the inverse matrix of g. 

In the situation at hand, we simply choose the conformally flat Riemannian metric 
g = p~ 1 {dx) 2 on the Euclidean unit ball -£>i(0), and define the weight to = p a+N A. 


With these ingredients, (4.1) becomes the heat equation 


d t w - A^w = f. 

Notice that by setting s = arcsin(|x|), the metric transforms into 
^g = ^p _1 (da;) 2 = (ds) 2 + (tan s) 2 (d£ S N-i) 2 , 

with d£ S N~i being the length element on the unit sphere. This representation in¬ 
dicates that the manifold degenerates at its boundary where lim s _ ) .| tans = oo. 
Loosely speaking, g can be considered to be “half way” between the Euclidean met¬ 
ric ( dx ) 2 and the metric on the hyperbolic Poincare disk p~ 2 (dx) 2 . We expect a 
deeper understanding of M. from the study of geodesic curves conducted in Subsec¬ 
tion 14.11 below. 

Our program for this present section is the following: In Subsection |4.1[ we compute 
the geodesic distance on (JA , g) as a function of coordinates on L>i(0). In Subsection 
4.3[ we study the homogeneous heat equation and derive Gaussian estimates for the 
heat kernel. These estimates can be used to establish maximal regularity estimates 


for the heat semi-group. Finally, Subsection 4 A contains the maximal regularity 
estimates for the inhomogeneous equation. 


4.1 Intrinsic distance, balls and volumes 

In the following, we study the geodesic distance on the (weighted) manifold intro¬ 
duced in the introduction of the present section. It is convenient to express this 
distance in terms of the conformally flat coordinates, that is, we study the distance 
induced by the metric g = p~ 1 (dx) 2 on the Euclidean ball Ri(0). In view of the 


interpretation of (4.1) as a heat equation on (A4,g,/i), this distance function on 


Ri(0) can be thought of as the intrinsic distance for diffusions of the form (4.1), 
in the sense that its second power measures the typical time scale at which heat is 
exchanged between two points. 

The intrinsic distance for (4.1) is (modulo a factor of \/A) defined as the quantity 


d(zi, z 2 ) = inf (L(T) : T joins Z\ and z 2 } , 
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where 


if T(a) = z i and T(6) = 
equivalent to minimizing 


m = 


f J>L 

Ja 


dr , 


Z' 2 - It is a standard observation that minimizing L is 


£(r) = 


f b r(r)l 2 

Ja l-|r(r)P 


dr , 


because L(T) 2 < (b — a)E(T) by Holder’s inequality with equality precisely if T 

ir'|2 

is geodesic, that is YVjpp = const. The functional E is strictly convex and admits 
hence a unique minimizer. Parametrizing the minimizer T by arclength, the geodesic 
equation reads 

|r'| 2 = 1 - |T| 2 , (4.2) 

and the Enler-Lagrange equations are 




(4.3) 


It is not difficult to see that geodesic curves through z\ and Z 2 are confined to the 
two-dimensional plane that is spanned by the points Zi, Z 2 and 0 in R v . Upon a 
rotation, we may thus write T = (x 1 y, 0,..., 0)^ 6 RxRx R^" 2 , so that ( |4.2[ ) and 
(4.3) translate into 


(x'Y + (y'Y = 1 -x 2 -y 2 


and 


x 


x* 


y z 


x 


and 




y 


X* 


X- 


yZ 


These equations are solved by hypocycloids, i.e., traces of fixed points on small 
circles of radius r < 1/2 that roll along the interior boundary of the unit ball: 


x{t) 

y{t) 


(1 - r) cos ( J ) + r cos ( \j ——t ) , 


1 — r 


(1 - r) sin ( J ^—t ) - r sin ( \j ——t ) . 


1 — r 


Here we have rotated and, if necessary, flipped over the two-dimensional disk in such 
a way that the geodesics hit the boundary at (1,0) and (cos(27rr), sin(27rr)). 

In spite of the good understanding of geodesic curves, except from a few particular 
cases, it is difficult to explicitly calculate the intrinsic distance between two points 
on M. On the positive side, the geodesic distance d on (A4,g) is equivalent to the 
semimetric 


d{zi,z 2 ) 


_Hi ~ Z2I_ 

a/ p{z l) + \J p(z 2) + y/\z\ — Z 2 


dehned on Bi(0 ): 
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Proposition 4.1. It holds 

d(z 1 ,z 2 ) ~ d(z 1 ,z 2 ) 

for all Zi, z 2 E B 1 (0). 


In the remainder of the present section, we will mostly choose d over d. It is thus 
convenient to speak of d in the following as a distance (or even the intrinsic distance), 
even though it lacks a proper triangle inequality. 

In our proof of Proposition 4.1, we follow [ 36 , Chapter 4.3]. 


Proof. We start with a study of two special cases. For two points on a straight 
line through the origin, we chose r = 1/2, so that x(t) = cos (t) and y(t) = 0. For 
convenience, we suppose that Z\ = az 2 for some a E [0,1]. Then 


d{zi,z 2 ) = |arccos(|zi|) - arccos(|z 2 |)| ~ \ VP( z i) ~ y/p(z 2 )\. (4.4) 


It remains to observe that the latter expression is equivalent to d(zi, z 2 ) if z\ = az 2 . 
Now consider two points on the boundary, \zf = \z 2 \ = 1, say Z\ = (1,0) and 
z 2 = (cos(27rr), sin(27rr)). If R(t) denotes the distance from the origin parametrized 


by arc length (cf. (4.2)), then 


R(t) 2 = 1 — 4r(l — r) sin 2 


2^(T 


and thus d(z \, z 2 ) = 27 ry / r(l — r). On the other hand, a direct calculation yields 
\zi — z 2 \ = 2sin(7rr). Solving both identities for r then gives 

d(zi,z 2 ) ~ y/\z 1 -z 2 \. (4.5) 


In the general case, we fix two non-parallel, non-zero z±,z 2 E B i(0). We start with 
proving the upper bound 

d(zi,z 2 ) < d(zi,z 2 ). (4.6) 

We distinguish two cases. First, if ^J\z\ — z 2 \ < \Jp{zf) + ^p(z 2 ), we have 


Vp( z i) - Vp( z 2 ) 


< 2 


\p(z 1 ) - p(z 2 )\ 


\/ P( z i) I + \/ P( z 2) + yj\ z \ — z 2 \ 


< 


d(zi, z 2 ), 


where we have used the definition of p in the second inequality. We do note loose 
generality by assuming that p(z\) < p{z 2 ). Hence, by triangle inequality 


d(zi, z 2 ) < d ( z 1} pPz\ ) + d ( Z \, 
Hi / V Hi 


that is, we replace the geodesics by first moving from z\ along a straight line in 
direction of the origin until we reach the sphere of radius |1 - From there, we move 


along the hypocycloid towards z 2 . By (4.4) and the previous estimate, the length 
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of the first piece of this curve is controlled by d(z±, Z 2 ). The length of the curve on 
the hypocycloid is cruelly estimated against the Euclidean distance, that is 


d ( pp-1,-2 ) ~ 


M 

kW 1 - 


< 


Hi - z 2 


Vp( z 2) Vp( z 2) 


The latter is estimated by d(zi,z 2 ) because of the assumptions p(z±) < p(z 2 ) and 
y/\zi-z 2 \ < \/p{zi) + y/p{z 2 ). 

In the second case, if ^J~p(zl) + \Jp{z 2 ) < \J\ Z \ ~ ^ 2 1, we use the triangle inequality 
d(zi,z 2 ) <d(z 1 , T ~ _ r ^ + d + d (-^^ 2 ') , 


Hi 


-1 *2 


z 2 


that means, this time our competitor curve first runs straight to the boundary and 
then along the hypocycloid to connect the projection points ph and pp Hence, 

applying fl4~4] ) and fl4~5| , we estimate d(zi,z 2 ) < p p{zi) + p{z 2 ) + ppi - z 2 \. 

The latter term is controlled by d(zi,z 2 ) thanks to the above assumption. This 
concludes the proof of (4.6). 

We finally turn to the proof of the opposite estimate 

d(z 1 ,z 2 ) < d(zi, z 2 ). 

We suppose that p(z 1 ) < p(z 2 ). Our Erst goal is the lower bound 


(4.7) 


d(z \, z 2 ) > min 

h>\/pb 2) 


*P - Vp( z 1) + 


z 1 - z 2 


*P 


(4.8) 


Indeed, on the one hand, we have the trivial estimate 

Hi - z 2 \ 


d(z 1 ,z 2 ) 


> 


ip* 


where ip* = sup 4 a/ p(T(f)). On the other hand, via the fundamental theorem and 
because T(f) e £>i(0), 


SUp pi ~ |T(f)| 2 < y/l - \zi\ 2 + 


TO I 


Vi - |r(f )| 2 


dt, 


and thus ip* — a Jp{z\) < dpi, z 2 ). Since ip* > a Jp{z 2 ), we have thus proved (4.8). To 
deduce (4.7), we again distinguish to cases. If the optimal ip satisfies the estimate 
ip > pPpi) + ppfo) + y/\zi ~ z 2 j, then 


d{z\.z 2 ) < a/ pi - z 2 \ <ip - pppi) ~ \fp(z 2) ib dpi, z 2 ). 
If not, we conclude 


d(zi,z 2 ) < 


|Q ~ z 2 \ 

ip 


d(zi,z 2 ). 


This proves the lower bound (4.7). 
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The knowledge of geodesics and of the intrinsic distance improves our understanding 
of the Riemannian manifold A4 . The boundary is at a finite distance from the 
interior of the ball. The shortest distance is realized by straight lines and is exactly 
d(z,dB) = arccos(|z|) ~ \Jp(z). By contrast, the length of curves parallel to the 
boundary diverges. For instance, the length of a circle of radius R is fy==. 

The study of the intrinsic distance function reveals the diffusion time scales on AT 
On a general manifold, diffusion over a distance d happens in a time t of order t ~ d 2 . 
In our situation, geodesics hit the boundary orthogonally. Hence, diffusion near the 
boundary happens primarily in normal direction and distances scale d ~ \JdEuci 
with the Euclidean distance dsuci • The diffusion time scale close to the boundary 
is thus t ~ dEuci■ On the other hand, in the interior where d ~ dsuch the diffusion 
time is as in the Euclidean setting, t ~ d 2 Eucl . The role of different diffusion time 
scales with regard to the actual position relative to the boundary will be reflected in 


our Calderon-Zygmund theory developed in Subsections |4.3| and |4.4| below. Notice 
that close to the boundary, the diffusion part — pA and the drift part (a + 1 )z ■ V 
of our linear operator are comparable. Here, the fact that a + 1 > 0 turns out to 
be of importance: Mass is transported towards the boundary and thus boundary 
conditions are not required. On the contrary, we will see that solutions satisfy 
some natural boundary conditions, which were described as “asymptotic boundary 
conditions’" 


m 


See also Remarks 4.1 and 4.2 in Subsection 4.2 below. 


The relevance of intrinsic distances for the regularity theory of the (linearized) 
porous medium equation was first exploited by Daskalopulos and Hamilton na 
and Koch [36] . 

We conclude this subsection with a study of intrinsic balls and their volumes. An 
open intrinsic ball of radius r around z is defined by 


B r (z) := \ y e B x { 0) : d(z, y) < r 


The first statement of Lemma 4T below shows that such a ball is comparable to the 
Euclidean ball with the same center z but radius r (r + \Jp(z ) 


Lemma 4.1. There exists a constant C > 0 such that for all r > 0 and z G RfyO), 
it holds that 


B 


r+VpW)^ 0 C Br ^ C B Cr(r+^j)( Z ) R Bl ^' 


Proof. It is enough to study the case r < 1. We fix z,z' G RfyO). Suppose that 
d(z, zJ) < r. We will first establish the second inclusion by showing that 


A < 


+ 


(4.9) 


From the definition of d we notice that 

\z - z'\ < r (\/p(z) + \Zp(z') + \J\z - R|) , 
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which with the help of Young’s inequality turns into 


i* - a % r + vp ( z '))+ r2 


\z , \ 

If \z \ > |z|, then the claim follows. Otherwise, if \z\ < \z\, we set z Wz. 
We recall that the shortest curve connecting z and z lies on a straight line and 
this line segment realizes the minimal distance betw een z and the ball of radius 

d(z,z) < d(z,z') < r. On 


4.1 


\z\ around the origin. Hence, using Proposition 
the other hand, because the geodesic curve is given 
[arcsin(|Y|), arcsin(|z|)], it holds 


d(z,z ) = arcsin(|^|) — arcsind^')) = 


"h'l 


'N 


VT 


T- 


yy T(t) = A sin(i) where t G 


dr ~ ^p(z') - y/p(z). 


We thus deduce that \Jp(z') < \Jp(z) + r, and therefore (4.9) follows. 

To prove the first inclusion, we assume that \z — z'\ < r + a/ p{z) s j . Then, by the 

definition of d and monotonicity 


d(z , z) < 


r r 


+ \fp(p 


VKz) + Vp( z ') + J r ( r + y/p&j) 


< r. 


This concludes the proof of Lemma 4.1 


Towards the end of this subsection, we will gather some technical results that prove 
to be helpful in the subsequent analysis. The first result provides two comparison 
formulas for balls whose centers are either relatively close or relatively far away from 
the boundary. 


Lemma 4.2. There exists C > 0 such that for all z G -£>i(0) and r > 0 the following 
holds: 

1. If z 7 ^ 0 and a/ p(z) < r, then 


B r 2 — nfii(O) C Bq t (z) and B d (z) C B c r 2 t~t nBi(O), 


and a Jp{z') < r for all z' G B d (z). 
2. If a Jp{z) > r, then 


B 


\/pb) 


(z) n ^(0) C B d {z) and B d {z ) C B 


CryfpfA) 


»nB 1 (o), 


and p[z !) < p(z) for all z! G B d (z). Moreover, if a/ p(z) > 2 Cr, then also 
p(z') > p(z) for all z' G Bf(z). 
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Proof. The inclusions in the first statement are immediate consequ ence s of the trian¬ 
gle inequality for d, the fact that d(z , z/\z\) ~ \Jp(z), and Lemma 
for any z' G Bf:(z), it holds 


4.1 


For instance, 


d (j, < d{z', z) + d (^z, < r + y/pC z ) < r, 

which implies the second inclusion via Lemma 4.1 The first inclusion can be estab¬ 
lished analogously. Moreover, ^fp < r in Bffz) follows from the second inclusion. 
The second statement is an application of Lemma 4.1 and the assumption. Indeed, 
the comparison formula for the balls follows from Lemma 4.1 by further estimating 
the radii. Then, by triangle inequality we have for all z' G B^(z) C B Cr ^^—(z) 

that \z\ < Cp(z ) + \z'\ and thus p(z') < p(z). Moreover, using \z'\ < Cr^/p(z) + \z\, 
we obtain 

1 - \z'\ >l-\z\- Cr\fp(z) > ^p(z) + ^y/p(z) - 2 Crj . 


For \/ r p(z) > 2 Cr this implies p(z') > p(z) as desired. ■ 

The next lemma provides estimates on the intrinsic volumes of balls, where “in¬ 
trinsic” refers to the volume form inherited from the weighted manifold discussed 
in the introduction to this section. To fix a notation reminiscent of a, we set 
p a := p u C N L -£>i(0) (= co volg), where £ iV denotes the Lebesgue measure, and for 
any measurable A C R jV , we write \A\ a = p a (A). Notice that for any a > — 1, 
p>(j defines an absolutely continuous finite Radon measures, which has precisely the 
same null sets as C N L £>i(0). 

Lemma 4.3. For any r < 1 and z G -Bi(O), it holds 

| B?(z)\ a ~r N (r + ^fa)) N+2ff . 


Proof. The statement is an immediate consequence of Lemma 4T 
We finally have: 


Lemma 4.4. Let r > 0 and z, zJ G Ri(0). Then it holds 




< i + 


d(z, z'] 


max{A^,2AT+2(j} 


In particular, 

. max{ AT,2iV-|-2cr} 

r + VpR) < L , d(z,z') \ 
r + y/p(z') ~ v r ) 

Proof. The results are trivial if r > 1. Otherwise, if r <C 1, we invoke the triangle 
inequality to verify the inclusion B^(z) C B^ r+d , Z ,^(z'). An application of Lemma 
El then yields the desired estimates. ■ 
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4.2 Existence, uniqueness and energy estimates 


The objective of this subsection is to give a precise meaning of a solution to (4.1) 


and the corresponding initial value problem, state existence and uniqueness results 
and derive first regularity estimates in the natural Hilbert space setting. 

We start with some considerations from a functional analysis point of view. The 


(sub-)elliptic operator occurring in (4.1), C := —pA + (a + 1 )z ■ V, is nonnegative 
symmetric with respect to the L 2 := L 2 (/i CT ) inner product, because 


(pCifi dp a = / V<p ■ V^dpcr+i for all <p,-0 G C' oo (I?i(0)). 


(4,10) 


As the space of test functions, C°°(i?i(0)), is densely contained in the weighted 
Sobolev spaces a+k , where 


m = Li 


H. 




= l: 


n H l a+ 1 n 


nH 




and H* +e = H e (p a+ £), cf. J4T, Lemma 2], C extends to a self-adjoint operator on 
L 2 (Friedrich’s extension). The following notion of weak solutions seems thus to be 
natural: 


Definition 4.1. Let 0 < T < oo and f G -/^((CfT);!! 2 ) and g G L 2 . We call w a 


weak solution to (4.1) with initial datum g, if 


wd t ( dfi a dt + 


Vw ■ VC d/j a+ idt = / g(dfj> a + 


ffd[i a dt (4.11) 


for all f G C°°([0, T ) x ^(0)) with spt ( C [0, T) x -BdO). 


Here and in the following, we use the convention that if domains of integration in 
the spatial variables are not specified, then we integrate over Hi(0). 


Remark 4.1 (Asymptotic boundary condition). Notice that the definition of weak 
solutions (and the one of C above) contains some (weak) information about the 
behavior of solutions near the boundary. Indeed, if w is a smooth solution, then 
necessarily 

lim p(zY +1 z ■ Ww(z) = 0. 

Mti 


This condition rules out functions that grow rapidly at the boundary. It is proved 
in the course of the present section (and announced in Theorem f.l above), that for 
sufficiently regular data, solutions are in fact C 1 up to the boundary. Such solutions 
trivially satisfy the asymptotic boundary conditions above. In fact, these boundary 
conditions are also the natural boundary conditions in the sense of calculus of vari¬ 
ations as (4.10) can be derived as the Euler-Lagrange equations for the (weighted) 
Dirichlet energy, cf. m Lemma 5], See also RemarkvfJf l 


Our first result concerns the wcll-posedness of the initial value problem in the Hilbert 
space setting. 
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Lemma 4.5. Let 0 < T < oo and f G L 1 ((0,T);L^) 
unique weak solution w to (4.1) with initial datum g. 
with w(0) = g and the energy identity 


and g G L 2 . There exists a 
Moreover, w G (7([0, T\]L 2 a ) 


1 

2 


«'Cn||£,+ 




fw dp a dt 


holds. 


Proof. The proof of such a result follows a standard procedure. For instance, ex¬ 
istence can be obtained via an implicit time-discretization and the elliptic theory 
derived in the appendix of m ■ The energy identity follows from choosing w as a 
test function in the weak formulation, which is admissible by the density of smooth 
functions in Hf a+l . Uniqueness is a consequence of the energy estimate and the 
linearity of the equation. ■ 


Remark 4.2 (No boundary conditions). It is remarkable that the weak notion of 


solutions in Definition guarantees uniqueness without specifying the boundary 
condition on <TE>i(0). This observation goes far beyoiid Remark f.l: In fact, the 
parabolic equation (4.1) is well-posed in L 2 only if boundary conditions are not 


prescribed. Roughly speaking, this phenomenon stems from the fact that close to the 
boundary, pA and z ■ V are both of the same order and the latter term drives mass 
towards the boundary. It seems that phenomena of this type have been systemically 
studied only very recently — at least in the unweighted setting, see 


We prove maximal regularity in L 2 for the equation with zero initial data. 


Lemma 4.6. Let 0 < T < oo and f G L 2 ((0, T); L 2 a (B)). Suppose that w is the 


weak solution to (4.1) with initial value g = 0. Then the mapping t K > || Vu^f)!^ 


is continuous in [0, T] and Vw(0) = 0. Moreover, 


<?+1 


/ \\d t w\\ 2 L2 dt + / \\Vw \\ 2 L 2 dt + 

1 0 Jo 


\V 2 w \\ 2 L 2 dt< 

(7 Z 


hdt. 


Proof. Let 0 < t\ < ^ < T be arbitrarily given. Using a regularized version of 
C — X[ti,t 2 \9tW as test function in (4.11), we infer the identity 


f {d t w) 2 dg a dt+^\\Vw(t 2 )\\ 2 L 2 +i 





fd t w dp a dt, 


from which we deduce that \\d t w\\ L 2 ^ L 2 ^ < ||/||z. 2 (z.g.), because g — 0 and Vw(0) = 
0. Moreover, further examination of the above identity reveals that [0, T) 3 t e->■ 
II Vu; (*)IIl * +1 e R is a continuous mapping. We deduce hence that it is enough to 
study the elliptic problem 


- p _<T V • (p° +1 Viv) =/:=/- d t w G L 2 ((0, T); L 2 ) (4.12) 
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pointwise in time, that is, we suppose that / € L 2 . Then the statement of the 
lemma is a consequence of the maximal regularity estimate 


\V 2 w\ 


Ll+2 


< 


Ll 


(4.13) 


and 


HViolla 


Ll- 


(4,14) 


Our arguments for (4.13) and (4.14) will be quite formal and can be made rigorous 


using standard approximation techniques and the density of smooth functions in 
Haa+ia+ 2 - We start by proving the auxiliary estimate 


IIVHI 


^ + i 


< 


Ll- 


(4.15) 


For this purpose, we test (4.12) with p a w and obtain using the Cauchy-Schwarz 


inequality and f f dp a = 0 (notice that spatially constant functions are admissible 
testfunctions) that 


\Ww\ 2 dp, a+l < f 2 dfi a (w — c) 2 dp c 


1/2 


where c is an arbitrary constant. Combined with the Hardy-Poincare inequality 


inf (vj — c) 2 du 0 

ceR / 


< 


\Vw\ 2 dp a+ i, 


which is proved in m Lemma 3], this implies (|4.15[). 


For (4.14), we first let r] be a smooth cut-off function that is supported outside B i (0) 
such that rj = 1 outside B i(0). Testing (4.12) with p a i ]^ • Vm then yields after 
multiple integration by parts 


^ ^ 1 / r)\Ww\ 2 dp a < [ rjf-^ ■ Vwdp a + C / | Vw\ 2 dp a+1 , 


and via Young’s inequality 


IBi(0)\Bi ( 0 ) 


r r (J4TT5J r 

\^w\ 2 dp a < / / 2 d/i CT + / | Vic| 2 dp a+ i “ / f 2 dp, 


On the other hand, because p ~ 1 in 5i(0), we get from (4.15) that 


'B 1 (0) 
2 


| Vin| 2 dp a < / f 2 dp 0 


The last two estimates together yield the desired bound (4.14). 
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We finally turn to the (formal) proof of (4.13), which is based on a suitable choice 


of testfunction and a multiple integration by parts. Indeed, testing with — p a+1 Aw, 
we derive 


V 2 w | 2 dfi a+2 


= — fAwdp a+1 + (a + 2) z ■ (V 2 w'Vw) dp a+l — / A wz ■ X7iudp a+ i, 


equalities and (4.14). 


from which we deduce (4.13) with the help of the Cauchy-Schwarz and Young in¬ 


We conclude this subsection with a symmetry observation for the homogeneous 
problem, that is / = 0. 

Lemma 4.7. If W\ and W 2 are two solutions to the linear homogeneous equation 


(that is, (4.1) with f = 0) with initial data g\ and g 2 , respectively, then 


w\(t)g 2 djicr — / giw 2 (t)dfjL 0 


for all t > 0. 


This lemma shows in particular that the semigroup operator e tc is self-adjoint on 

LI 

Proof. We define the time inversion operator S : [0, t\ —> [0, t] by «5(r) = t — r. Then 

J w\(t)g 2 dn a — J giw 2 (t) d/a a 

= [ w 1 (t)(w 2 oS)(t)dfj, cr - [ w 1 (0)(w 2 oS)(0)dfi a 



(d t wi(w 2 o S) - wi(d t w 2 ) o S ) dp a dt. 


Using the equations for w\ and w 2 and integrating by parts yields the desired result. 


4.3 Carleson measure estimates for the homogeneous prob¬ 
lem 

In this subsection, we derive regularity estimates for the homogeneous equation 

d t w - p- a V ■ {p a+1 Ww) = 0, (4.16) 

w( 0, •) = g. (4.17) 

We suppose that the initial datum is a Lipschitz function with respect to the Eu¬ 
clidean topology, that means 11£/11Lip = ||Vg|| £00 <C OO. 

In a first step, we derive local Lf regularity estimates. 
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Lemma 4.8. Let 0 < £ < i < 1 and 0 < <5 < 5 < 1. Let w be a solution to the 


homogeneous equation (4.16). If 5 and 5/5 are sufficiently small, then the following 
holds for any k £ N 0 , /3 £ and 0 < r < a/2: 

1. If z 0 £ B\ (0) is such that ■>/ p{zq) < 5r and r > 0, then 


rr-\-r z 


' r-\-er 2 


\\ d t^M\l2 {B f {zo ))dt<r 41/31 4k 

or 


r»T+r 2 


t+£T 2 


\ W \\Ll(Bl(z 0 )) dt - 


2. If z 0 £ B\ (0) is such that ^J p(z 0 ) > 5r and t > 0, then 


f*r+r z 


' T-\-er 2 


\\ d t d z W \\ll(B<l {z 0 )) dt ~ r ^ m p( Z o) 


r»r+r 2 


’ r+er 2 


\W\ 


L 2 ABf(z 0 )) 


dt. 


Of course, the constants in the inequalities (which are suppressed in our notation) 
depend on the particular values of k , /3, £, £, 5, and 5. For our subsequent argu¬ 
ments, only the dependence on r, r, and Zq is of relevance. We will thus neglect the 
dependency on the other variables for notational convenience. 


Proof. Upon a shift in the time variable, we may without loss of generality assume 
that r = 0. 

We first consider the case \j p(zq) < 5r. In view of Lemma 4.2 it suffices to prove the 
statement for z 0 £ dB i (0), provided that 5/5 is small enough. In the first statement 
considered here, we may thus replace the intrinsic balls by Euclidean balls with 
squared radii. 

We choose £ < i < £ and 5 < 5 < 5 and let g be a smooth cut-off function 
that is supported in (er 2 ,r 2 ] x Bp r 2 (zo) fl-Bi(O), and satisfies 77 = 1 in [ir 2 ,r 2 ] x 
Bjj 2 r 2 (z 0 ) D £i(0) and \d/dfrj\ < r~ 2k ~ 2 I /3 L For the arguments that follow it is con¬ 
venient to first consider the inhomogeneous equation 

d t w - p~ a V • (/U +1 VU7) = /, 


for some / that will be specified later. Then gw solves 
dtfgw) — p~ a 'V ■ (p cr+1 V(? 7 , u;)) = gf + {d t g — p~ a V ■ (p CT+1 V//)) w — 2 pS7g ■ Via. 


For further reference, we denote the left-hand side of this equation by /. Notice that 
g acts as a spatially constant function on the differential equation if 5 2 r 2 > 2. In 
this case, the equation and most of the arguments that follow simplify dramatically. 
We will thus focus on the case where g is not spatially constant. Testing with p a gw , 
using the properties of g and observing that p(z) < r 2 for z £ Bs 2 r 2 (z 0 ) D f?i(0), we 
derive the Caccioppoli estimate 



< 

r^j 



|Vwf dp a+ idt 


f 2, dprjdt + 


B s 2 r 2(z 0 ) 



W 


1 B s 2 r 2{zo) 


1 dp a dt. 


(4.18) 
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Here we have written Bi(zq) for Bg(zo) fl -Bi(O) and will do so in the following 
whenever it seems convenient. 

The estimate in (4.18) is not yet the statement we are looking for since the weight 
on the left-hand side is still too small. We choose a new cut-off function by replacing 
(e, S) by (e, 5) and (e, 5) by (e, <5) in the above definition. We denote this function 
again by p. Invoking the estimate from Lemma T6 with w and / replaced by rjw 
and / we obtain 


(( d t (pw )) 2 + \W(r)w)\ 2 + \p\/ 2 (pw)\ 2 ) dp a dt 


< 


j dp a dt. 


By the definition of /, using the properties of p and the fact that p(z) < r 2 in the 
support of p, 


f 2 dp a dt 


lir 2 JB- 5 2 r 2 {za) 


< / / (/ +~iW ] d’hadt + — 

Jir 2 J B/j 2 r -2 (zo) V r / r • 

Combining the previous two estimates now yields 

z* 7 * 2 r 

((d t w ) 2 + |Vrc | 2 + |pV 2 w| 2 ) dp a dt 


| Vw| 2 dp a+ idt. 


ler 2 JB t ,2 r 2bo) 


< 


Us 2 JB g 2 r 2(z 0 ) 


f 2 dp^dt + — 


r der 2 JB s 2 r 2 {z 0 ) 


w 2 dp a dt. 


(4.19) 


In particular, choosing / = 0, we obtain the statement in the cases (k, |/3|) = (0,1) 
and (k, |/3|) = (1, 0). 

To prove the analogous statement for the second-order spatial derivatives, we have 
to lower the exponent on weight factor. For this purpose, we first differentiate the 
equation in direction of the tangent vectors close to z 0 G dB. Without loss of 
generality we assume that z 0 = e\. We let R 2 ,..., Rn be rotation matrices in the 
sense that RjZ _L 0 for all i and all z G R A , and such that {z, R 2 z ,..., Rnz} forms a 
basis of R a for any z G Bt(e i), provided that £ is chosen sufficiently small. Suppose 
moreover that the R £s are such that R,ei = At any point z G B S 2 r 2 (z 0 ) we let 
di = Zi ■ V where z t = RiZ for some fixed i. Differentiating the parabolic equation 
yields 

dt&iiv - p _<T V ■ {p a+l Vd iW ) = -2 pRi : \7 2 w + (a + l)d iW , 


2 E.g., in 3 D, we let 


/ 0 

-1 

0 \ 


/ 0 

0 

- 1 \ 

1 

0 

- 1 

and R3 = 

° 

0 


V 0 

1 

0 / 


V 1 

1 

0 / 


Rn = 


so that R 2 e\ = e 2 and R^e i = e3. Then, for all 2 G R 3 it holds R\z _L z and R 2 z _L z. Moreover, 
for any z € BRe 1) with l sufficiently small, the set {2, R 2 z , -R32} forms a basis of R 3 . 
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where A : B := Ylij for any two matrices A and B. Applying (4.19) to this 

equation (with suitable changes between (e, h), (£, 5), and (£, 5)) and then once again 
to the original equation yields 


ler 2 J B- s 2 t2 {zo) 

< f 


(■ d t diW ) 2 + \'Vd i w\ 2 + \pS7 2 diw\ 2 ) dp a dt 


lir 2 JB- 5 2 r 2(z 0 ) 


\pS/~w\ dp a dt + — 


r 3er 2 JB- S 2 r2 (z 0 ) 


| Vw| 2 dp a dt 


x rr 2 


hr 2 JB 5 2 r 2(z 0 ) 


w 2 djlrjdt. 


Applying (4.19) again, we can change the order of differentiation on the left-hand 
side and obtain 



< 

r^j 


’ B s 2 r 2 (20) 
„G2 


did t w | 2 + \djVw\- + \pdi'V 2 w\ 2 
w 2 dp a dt 1 


ler 2 jB 6 2 r 2{z 0 ) 


dp a dt 


and the estimate is independent of the particular value of i 6 {2,...,N}. This 
procedure can be iterated and yields for any multi-index (3 e 



< 

r^> 


r i(\0\+i) I I 

Jer 2 J B s 2 r 2 (20) 


w 2 dp a dt. 


(4.20) 


Here, we have used the notation z' = (z 2 ,..., zy) T and defined <9f, in the obvious 
way. 

In the following, we will write V' to denote the tangential gradient operator, that 
means, (<9 2 ,..., d]\r) T . 

We turn to the estimates of the derivatives in direction normal to the boundary. 
Differentiating the equation in normal direction, using d\ — z ■ V yields 


dAw - p~ a ' V ■ (p a+1 Vd\w) + Aw = 0. 

Observing that z ■ V d\w = d\ w + (z — z) <E> (z + z) : V 2 w + z <E> z : V 2 u>, where 
5 = z/\z\, and A w — z® z : V 2 w = 0(|V'Vw|) because {z, Riz,... Rnz} is a basis 
of R a and R, z _L z, 0 we can rewrite this equation as 

dAw - p" (,T+1) V • {p a+2 Vd lW ) = &xw + O (|pV 2 w|) + O (| V'Vw|) . (4.21) 

3 Since {z, R 2 z ,..., Rnz} is a basis of H N , we find \j G R such that e, = \\Z + X)^=2 KjRj z - 
In particular, for any matrix A G R WxJV , we have that An = 6 i®ei : A = (X a z + Y 3 jL 2 A ijRjZ J <E> 
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As in the derivation of (4.19), we can smuggle a cut-off function into the equation. 
Applying (4.19) with / being the right-hand side of the resulting equation and with 
cr replaced by a + 1 yields the estimate 


'er 2 JB S '2 r 2(zo) 

< r 2 ( 


(d t diw) 2 + |V<9iu>| J + \pV 2 d 1 w \ 2 ) dp a+ idt 


'er 2 J B s 2 r 2 bo) 

i r 2 r 


\pX7 2 w\" + IV'Vwf) dp a+ idt 




Ur 2 JB s 2 r 2 (z 0 ) 


\pWdiw\ 2 dp, a+1 dt + — 


hr 2 JB 5 2 r 2(z 0 ) 


| Viu| 2 dp a+ idt. 


'er 2 JB- S 2 r 2 {z 0 ) 

< i r 2 

r 


The terms on the right can be estimated using that p a +i < min{l, r 2 }p, a in the 
domain of integration. We change the order of differentiation and obtain via (4.20) 

(<9i d t w) 2 + <9iV'«;| 2 + \pd{\7 2 w\ 2 j dp a+ idt 
i / w 2 dp a dt. 

er 2 JB s 2 r 2 {zo) 

We can argue similarly for the third order derivatives. Differentiating once more in 
transversal direction yields the estimate 


1 er 2 J B- S 2 r 2{z 0 ) 


(■ d'fwY + I df\7w\ 2 dp u+2 dt < — 


Ur 2 JB g 2 r 2(z 0 ) 


w 2 dpvdt. 


We finally use a Hardy-type inequality to control the second order derivatives with 
respect to the measure p a : Using the fact that (1 — \z\ 2 ) a ~ —V • (z(l — \z \ 2 ) a+1 ) 
close to the boundary, it can be proved via integration by parts that 


Ur 2 j B~ s 2 r 2{zo) 

< r 2 r 


' er 2 JB 5 2 r 2(z 0 ) 


did t w\ 2 + \di'Vw\ 2 j dp„dt 

\dld t w\ 2 + |<9 2 Vrc| J ) dp u+2 dt. 


+ YljL 2 ’■ A. If A is symmetric, an expansion yields 


N N 

A* i = (Ai i) 2 z <8> z : A + 'y ' A^A ik(Rjz) ® (Ruz) : A + 2 ^ ' XijXik(Rjz) 0 z : A. 

j,k —2 j,k =2 


Notice that a ■ z = Xu by orthogonality, and thus 1 = |z| 2 = Y2iLi( e i ' z ) 2 = Eili(^i) 2 - In 
particular, summing over all i yields 

N N N N N 

tr (A) = ^A ii = 5i8)0:A + ^ ^ XijX ik (Rjz) 0> (Rkz) : A + 2^ ^ X l jX ik {Rjz) <g> z : A. 

i—1 i =1 j,k—2 i—1 j,k=2 
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A combination of the previous estimates now yields 


\did t w\~ + \diViu \ 2 ) dp a dt < — 


hr 2 JB s 2 r 2{zo) 


w 2 d/J, a dt , 


and thus, together with (4.20) we get the first statement of the lemma for (k, |/3|) G 
{(1,1), (0, 2)}. Control on higher order estimates can be obtained via iteration. 

The second statement can be proved in a very similar (but easier) way. Indeed, 
choosing <5 smaller if necessary, Lemma T2 enables us to derive the statement in the 
Euclidean and unweighted setting. (In other words, the equation is strictly parabolic 
in this situation.) Such a result being standard, we omit further details. ■ 

In the following lemma, we establish uniform bounds on solutions and their deriva¬ 
tives by combining the estimates from the previous lemma with Morrey-type in¬ 
equalities. 


Lemma 4.9. Let w be a solution to the homogeneous equation (4.16) and let zq G B 
and r > 0. Let 0 < £ < 1 and 0 < 8 < 1 sufficiently small. Then, for all k G N 0 , 
/3 G and 0 < r < 1, the following holds: 


1. If a / p(zq) < dr, then 


I 2 < 


r ~4k-4\f)\ rr+r 2 


r 2 \B?(z 0 )\a 
for all ( t,z ) G [r + £r 2 ,r + r 2 ] x B$ r (z 0 ). 
2. If \/p{zq) > dr, then 

r -4k-2\p\ 




w 2 dp a dt 


(4.22) 


fT+r z 


\dfd^w(t,z)\ 2 < 


r 2 \B^(zo)\ cr 
for all ( t,z ) G [r + £r 2 ,r + r 2 ] x B$ r (z 0 ). 


iBd(zo) 


w 2 dp a dt 


(4.23) 


Proof. Upon a shift in the time variable it is enough to consider the case r = 0 
only. As in the proof of the previous lemma, we will often write the handier B^Zq) 
instead of B^zf) fl L?i(0). 

We will start considering the scenario p(zq) < dr. According to Lemma 4.2, it is 
enough to prove (4.22) in the situation where Zo is located on the boundary of the 
unit ball. In this case Bf r (z 0 ) can be replaced by B 5 2 t 2 [zq). The statement in (4.22) 
now follows from the previous lemma via the estimate 


1 N+p 


HvirsEE 


r U+4\^\ 


£=0 | 7 |=0 


r 2 \B 6 2 r 2 {zo)\c 


'er 2 j B s 2 r 2(zo)nB 


(df d^v) 2 dp a dt, 


setting v = dfdfw and where p G N 0 is such that 0<cr + l— p < 1 and (t, z) G 
[er 2 ,r 2 ] x B S 2 r 2 (z 0 ). (We do not attempt to give a sharp estimate here.) Since z 0 is 
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on the boundary, if 5 is sufficiently small, we may replace -Bi(O) by the half-space 
H = {z £ R N : Zn > 0} and assume that z o = 0. By the equivalence of norms, we 
may also consider cubes instead of balls, and upon rescaling t = r 2 t , z = r 2 z and 
v — v, the previous estimate follows from 

1 N+p x 

Hi z )\ 2 < / / ( d i d Jv ) 2 *N dzdt. 

£=0 |<y|=o d £ dQ s2 {0)nH 

Here, Q r (0) denotes a cube in R Af , centered at 0 and with side length 2 r. In view 
of the elementary inequality sup [ab] |/| < \\f\\ L Ha,b) + \\f'\\L^a,b) (with constants 
dependent on a and 6), we further reduce the statement we have to prove to 


\v[z N , 


max{p,l} g 

12 \ A ~ y \2 




(4.24) 


The latter is established as follows: We drop the hats and the subscript N from z. 
For z £ [0,5] and q = a + 1 — p £ [0,1), we have 

/ 5 p5 / pS p6 \ 1/2 / z 1 ! i \ 1/2 

|u| dz + J \v'\dz< (J v 2 z q dz + J ( v') 2 z q dz\ (J —dz\ , 

where we have used the Cauchy-Schwarz inequality. Since q < 1, the latter integral 
converges. Notice first that if p — 0, that is a < 0, we can lower the weight, 


z q = z a+1 < z a , and (4.24) follows from the previous estimate. The same holds true 


if p — 1, because then q = a. From now on we assume that p > 2. Letting rj denote 
a cut-off function supported in [0, 25} with rj — 1 in [0, 5], we obtain via integration 
by parts 

[ pf 2 z q dz< [ rj(f) 2 z q+1 dz+ [ (rj + \rj'\)(f) 2 z q+1 dz. 

Jo Jo Jo 

Choosing f = v and / = v', and repeating this procedure p — 1 times yields 


(v 2 + [y') 2 )z q dz < 


v 


-25 


(d™v) 2 z q+p ~ l dz. 


By the choice of p and the definition of q, it is q + p — 1 = cr, which proves (4.24). 


We finally comment on the proof of the second assertion. Thanks to Lemma 4.2 


the assumption p(zo) > 5r converts the problem into the unweighted Euclidean 
setting. Thus, a standard Morrey inequality is applicable and yields the desired 
result. We omit details. ■ 


For two arbitrarily fixed real numbers a and b, b > 0, we define the auxiliary function 
Xa,b ■ B 1 (0) x -E>i(0) [0, oo) by 


Xa,b(z, Zq) 


ad(z, z 0 ) 2 
\Jb + d(z, Zq ) 2 
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where 



\z - Zq\ 


p( z ) + p( z o) + 

\z - Zq\ 


d(z, z 0 ) 2 := . . -L - r ~ d(z,z 0 ) 2 . 


A short computation shows that 

v75)|Ww(*,*l)| £ |o|, 


(4.25) 


which in particular implies that Xa,b(',z< o) is Lipschitz continuous with respect to 
the induced metric d, that is 


\Xa,b{z,Zo) ~ Xa,b(z', z 0 )\ < \a\d(z,z'). 


(4.26) 


Indeed, if we denote the gradient on (A4,g) by V g , then Lipschitz continuity of 
a mapping y : A4 —)■ R with Lipschitz constant proportional to |a| follows from 
the gradient estimate a/ g(V g y, V g y) < |a|. Since in our case g ~ p~ 1 (dz) 2 and 


V g ~ pV, the latter is equivalent to (4.25). 


The function Xa,b will be a key ingredient in the derivation of the Carleson measure 
estimate (Proposition |4.2[ ) and the Gaussian estimate (Proposition |4.3[ ) below. The 
following result is an intermediate step towards both estimates. 


Lemma 4.10. Letw be the solution to the initial value problem (4.16), (4.17). Then 
there exists a constant C > 0 such that for all z, zo G £>i(0), 0 < r < a/2 and t > r 2 
and any k G N 0 and (3 G N/ it holds 




’0 

2fc—I 


(r + y/rfz))- 


\B*(z)\ 


1/2 


_ e Ca 2 t-Xa,6G^°)|| e Xa,(,(- ,2°)^|| l2 


Proof. For abbreviation, we will write y := y a ,b( ■ ,Zo). Notice that e x w solves the 
equation 

d t (e x w) - p~ a W ■ (p a+1 W(e x w)) = -p~ a W ■ (p CT+1 Ve x ) w - 2 p\7w ■ Ve x . 

Hence, testing with p a e x w and integrating by parts a multiple times yields the 
identity 

“ [(e x u>) 2 dp, a + I \V(e x w)\ 2 dp a+1 = f | Ve x \ 2 w 2 dp a+1 . 


Invoking the gradient estimate (4.25) and a Gronwall argument, this implies 

f (e x w{t)) 2 dp a < e 2Ga2t [(e x g) 2 dp a 


(4.27) 


for some C > 0. We now combine this estimate with the L°°-L 2 bounds from 


Lemma 4.9 Estimates (4.22) and (4.23) together imply the bound 


-4fc-2|/J|/ , /T7 Wa-2|/3| 

m w ^< r - — 1 ( ;+/; ()) 


sup e 2x< -^ sup 


z&B${z) 


t<T+r 2 JB^{z) 


(e x w(t)) 2 dp 0 


39 



















for all 0 < r < 1 and t G [r + £r 2 , r + r 2 ]. Hence, applying (4.27), we have 


\d k t d 0 z w{t,z)\< r ^ + 1/31 -- x( " )+Ca2(T+r2) H^ 


mz)\ 


snp e 

zS-B^O) 


eX g\\ Li¬ 


lt remains to invoke the Lipschitz bound (4.26) and we choose t = r + r 2 for 
conclusion. ■ 

The above estimate is valuable for estimating the small time behavior of solutions. 
For large times, we can do much better by taking into account the spectral infor¬ 
mation on the linear operator obtained by the author in B33- We will recall the 
results of m in every detail in Section [6] below. At this stage, only the following is 
relevant: A solution w to the initial value problem (4.16), (4.17) can be written as 
the semi-flow w(t) = e~ m g on the homogeneous Sobolev space Hf +1 . Here TL is the 
self-adjoint (Friedrichs) extension of —pA + (cr + l)z ■ V to LT/ +1 . The regularity 


obtained in Lemma 4.10 guarantees that solutions are in the domain of PL for pos¬ 
itive times. The work Bn contains a computation of the purely discrete spectrum 
of PL which lies in the set (0, oo). In particular, if we denote by Ai the smallest 
eigenvalue, we have the spectral gap estimate 


Ww-PJPLw d/i a+ 1 > Ai / | Vu>| 2 d/i a+ i. 


(4.28) 


Notice that the latter is not true for the Hilbert space L 2 . Indeed, the latter contains 
nonzero constants (which are modded out in the diagonalization [47]), and thus, the 
corresponding operator on L 2 has, in addition, a zero eigenvalue. 

With these preparations, we are now in the position to state and prove the following: 


Lemma 4.11. Let w be the solution to the initial value problem (4.16), (4.17). Let 


k G N 0 and (3 G N(/ be such that k + \/3\ > 1. Then, for all t > ^ and z G -£>i(0) it 


holds 


|af8Mi,z)l<e- A “||V 9 |U 


£T + 1 


Proof. Since in is a solution to the homogeneous equation d t w + TLw = 0, it holds 
that 

dt2 / l Vu, | 2 ^+i = “ / ' VHwd/v+i. 


The spectral gap estimate (4.28) and a Gronwall argument then yield 

-- At l|Vslb 


||Vw|| L 2 +i < e 


K+i- 


Notice that c = f w dp a is a conserved quantity for the homogeneous equation 


(4.16). Hence, with the help of Lemma 4.9 in which we choose r = l ,£ = 1/2 and 
z o = 0, and a Hardy-Poincare estimate (e.g., m Lemma 3]) we obtain for every 
(t,z) G [r + 1/2, r + 1] x i?i(0) that 


\a^w(t,z )\ 2 < 


f*T +1 


r»T+l 


(w — c) 2 dn 0 


< 


\Vw\ 2 dn a+ i. 


We easily infer the statement of the lemma. 
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The first main result of this subsection is the control of the Carleson measures 
||w||x( P ) hi terms of the Lipschitz norm of the initial datum. 


Proposition 4.2. Let w be the solution to the initial value problem (4.16), (4.17). 
Then 

IM|A'( P ) + IMlLip < 115 1 11 Lip • 

Proof. Upon a shift in the time variable, the exponential decay estimates of Lemma 


4.11 imply that 


\dtd^w(t,z)\ < e Xlt \\w(l/2) || Lip , 

for ant t > 1, and thus the control of large time terms (r > 1) in the Carleson 
measures follows from the control of the small time terms (r < 1). 

Let z 0 G -Bi(O), 0 < r < 1 and t G ^y,r 2 ^. Observing that w — g(z 0 ) is a 
solution to the homogeneous equation with initial datum g — g(zo) and using that 


Xa,b( z o,Zo ) — 0, we deduce from Lemma 4.10 that 

MU, (w(t,z)-g(z 0 ))\ 

< r- 2l, - |gl (r+y/^))-W CaH [,.„i| 

m~„)P 


-^o|||l 2 11S'11Lip 


for all k G N 0 and f3 G N^. We now specify the choice of a and b by setting a = — 4, 
b = r 2 . We claim that 


I X_1 _2 (•>*<>) 1 
e r’ r 


-ZoWIl* < r(r + Vp(z 0 ))\B^z 0 )fJf 


(4.29) 


Prior to establishing this estimate, we derive the statement of proposition. Thanks 


to (4.29), the estimate from Lemma 4.10 (now with the above choice of a and b ) 

|3?afto(t,*>)| <r 1 - 2k - m (r+ \7W) 1 - | ' s| || 9 ||Lip, 


becomes 


for all £ 0 G -£>i(0), provided that k + \f3\ > 1. In the case (k, |/3|) = (0,1), the terms 
depending on r drop out and we obtain uniform control on Vw and thus on the 
Lipschitz norm of w. Otherwise, if (k, |/3|) G {(1,1), (0,2), (0,3)}, integration over 
Qr(z) and rearranging terms yields the desired estimates on d t Ww, V 2 w, and V 3 w. 
Notice that in this argument we additionally use the estimate r < r + ■>/ p(z o) < 


4.2 


a/ p(z) for all z 0 G Bf(z), which is a consequence of Lemma 
We finally turn to the proof of (4.29). For j G N consider the annuli Aj := 
Bj r (z 0 ) \ Bfj^yfeo), which cover i?i(0) and satisfy Aj = 0 if j }. An ele¬ 
mentary computation shows that X-i,r 2 ( z i z o) — ~^ f° r z ^ Aj. Hence, for all 


i < i. 



i{z,zo) 


z - z 0 \ 2 dp a (z) 


< 


e vPfr 2 


(jr + V P(z 0 )) 2 \Aj 
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by the comparison formula for Euclidean and intrinsic balls, cf. Lemma |4dj Invoking 
the volume formula for intrinsic balls from Lemma 14.31 we further estimate 


I< (jr) N (jr + Vp( z o)) N+2a ;$ J K \Br(zo)\a 
for some k > 0. Notice that the latter is only valid because j < -. Hence, combining 


the previous two estimates and summing over all j £ N 0 yields (4.29) as desired. 
This concludes the proof. ■ 


From Lemma 4.10 , we also deduce the existence of a heat kernel, that is, the integral 
kernel of the heat semi-group operator e~ tc (or e tA, ‘, if the stick to the notation of 
the introduction to this section). It is well-known that heat kernels are smooth and 
symmetric functions and inherit the semi-group property from the semi-group oper¬ 
ator, cf. 1261 Theorem 3.3]. For the sake of a self-contained presentation, we derive 
the properties that are relevant for subsequent analysis directly from the Lemmas 


4.7 4.10 , and 4.11 above. In particular, we establish sharp decay estimates on the 
kernel, that are called Gaussian estimates in analogy to the Euclidean situation. For 
general parabolic equation, such estimates we derived by Aronson in [5j. 


Proposition 4.3. There exists a function G : (0, oo) x L?i(0) x L?i(0) — y R with the 
following properties: 


1. For any g £ L 2 , the solution w to the initial value problem (4.16), (4.17) has 
the representation 

d£d%w(t,z)= [ d^d z G(t, z, z')g(z') dp a (z') 


for all k £ N 0 , (3 £ , t > 0 and z £ L>i(0) 


2. For every (■ t,z,z') £ (0, oo) x Hi(0) x Hx(0), it holds 

G(t, z, z 1 ) = G{t 1 z !, z), 
i.e., G(t, •, •) is symmetric for every t > 0. 

3. It holds 

d t G-p-°V z -(p° +1 V z G)= 0 


and, for all z £ Hx(0) 

p a G(t , z, ■) —* S z as t fO in the sense of distributions. 


In particular, G is a fundamental solution of (4.16). 
f. For all z,z' £ -Bi(O), t £ (0,oo) and k £ N 0 , f3 £ N A , it holds that 


\dtdgG(t, z, z')\ < 


Vi 


-2k-\p\ 


Vi + V~pV) 


\B d Jz)\l /2 \B d Jz')\l /2 


3 — Cd(z,z') 2 /t 
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5. For all z, z! £ -Bi(O), t £ (1, oo) and k £ N 0 , (3 £ with k + \/3\ > 1 it holds 

\d^G(t,z,z')\<e- x “. 


Proof. For the existence of the heat kernel, we observe that Lemma |4.10| applied 
with a = 0, and thus x = 0, shows the boundedness of the linear mapping 

Ll3 w(t,z) £ R, 


for all (t,z) £ (0, oo) x Bi(0). Hence Riesz’ representation theorem yields the 
existence of a unique function G k> p(t,z, ■) £ L 2 such that 

dfdgw(t, z) = j G kj p(t, z, z')g(z') dg a (z'). 

By uniqueness, it is clear that d^dfG = G k ,p, where G = Go,o- Moreover, symmetry 


follows from the symmetry of the heat semi-group operator in Lemma 4.7 


We now prove the Gaussian estimate. It is a consequence of the auxiliary estimate 

r -2k-\p\/ r , ./F(7))-\P\ 

I dtd?G(t, z, z')\ < - ,,T — - „ e C ° ■»>, (4.30) 


which holds for any z,z',z 0 £ L>i(0), t > r 2 /2 > 0 and a £ R, b > 0. Indeed, 
choosing zl = z 0 , b = 3 d(z, z') 2 and a = —2£ with i > 0, we deduce from (4.30) that 


|a‘0fG(M,X)| 




< 


\r + \J\ |z|) ^ Ace 2 t-id(z,z') 


mz)\i /2 mz')\ 


1/2 


Optimizing the exponent with respect to i yields it ~ d(z,z') ~ d(z,z'). The 
Gaussian estimate follows by choosing r = y r t. 

We turn to the proof of (4.30). Keeping a £ R and b > 0 be arbitrarily fixed, we 
write x — Xa,b( • , £o) for abbreviation. Let f £ Lf \= L 1 (/x (T ) be fixed and such that 


.9/ := e x \B r 


, 1/2 


£ £ L 2 . Then given by 


wt(t,z) = j G(t, z, z')g^(z') dn a (z') 

is a solution to the linear homogeneous equation with initial datum g^. Applying 
Lemma 4.10 with a replaced by —a, i.e., x replaced by — x, and with initial time 
t/ 2, we obtain 


\B r G)\a 


e x vj£ 


\Ll- 


(4.31) 


In particular, since g £ L 2 precisely if h — e x g £ L 2 , there exists a well-defined 
operator 


A : L 2 a 3 h (->• e x h = g t-)- w i —> w ( -, • ) i —> \B 


1 1/2 e x w[-,- ) £L c 
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This operator is linear and bounded. Indeed, thanks to Lemma 4.10 it holds 


or, equivalently, 


|B r i ( Z )iyVW|« l (|,d|<e c “ 2 ‘ /2 


\\Ah\\ L ~<e c ^ 2 \\h\\Ll- 


Ll-, 


By the definition of Ah and the symmetry of the heat semi-group operator, Lemma 
ED we have the identity 


J Ahh, d/J-cr = I e X W£ ■ J hd/j, a . 

Therefore, the action of the dual operator A* : ( L °°)* —>• L 2 a on functions ( 6 
with g^ G L 2 a is *4*£ = e~ x w% (|, •). In particular, because ||.4.|| = ||„4*||, we must 
have 

‘ A 'k < <k 2,/2 ||£|U>. 


\e x vj£ 


We now conclude 

J dt d z G (hz, ■)g^dn 0 

= dtd% Wt(t,z) 
T31 
W 


(4.32) 


\ 4 -A r 2k W(r + \Jp(z))~ 


\B${z)\ 


1/2 


e CaH/2+ X {z)^ e -x w ^ 


I Ll 


{05} r 2 k l^l(r + p{z))~ 


mz)\ 


1/2 


0 Ca?tl2+x(z) 


neik- 


The last expression in the above chain of inequalities is bounded for any £ G L* 
and thus, by approximation, the estimate holds for any such £. On the other hand, 
thanks to the duality L°° = (L^.)*, we have that 

| B r V) k 2 1 d,‘df G(M,z') | 

< sup //’e*|B/(-)|y 2 9 t ‘afG((,z, -)(,dn c , 


which in combination with the previous estimate yields (4.30). 

Since G is a smooth kernel, it satisfies (4.16) itself, and p a G(t,z, ■) —;> 8 Z in the 
sense of distributions. 

Finally, since G( ■, ■ ,z') is a solution to (4.16), an application of Lemma 4.11 and 
the Gaussian estimate yields 


|9?afG(MV)l <e- Al *||VG(l/2, -,z') |U; +1 < e-'“. 
This concludes the proof of Proposition |4.3[ 
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4.4 Carleson measure estimates for the inhomogeneous prob¬ 
lem 

This subsection is devoted to the study of the initial value problem for the inhomo¬ 
geneous parabolic equation 


dtw - p~ a \7 ■ (p a+1 Vw) = /, 

iw(0, ■) = 0. 


(4.33) 

(4.34) 


Thanks to the existence of a smooth solution kernel G(t, z, z') for the homogeneous 


initial value problem, cf. Proposition 4.3 above, solutions to (4.33), (4.34) have an 


integral representation. This is the content of the following 


Lemma 4.12 (Duhamel’s principle). For (t,t',z,z') G (0, oo) x (0, oo) x B i(0) x 
i?i(0) with t! < t define G(t, t'z, z') := G{t — t !, z, z'). If f G L 2 ((0, oo); L 2 ) and w 


is the solution to the initial value problem (4.33), (4.34), then 


w(t, z) = 



G(t, t\ z , z')f{tz') dp a (z')dt r , 


for all (t,z) G (0,oo) x i?i(0). 

Proof. The statement follows from the fact that G is a fundamental solution to 


(4.16), see Proposition 4.3 


In Lemma 4.6, we established L 2 (L 2 ) estimates on d t w, \7w and pS7w for solutions 


to the inhomogeneous problem ( |4.33| ), (4.34). As a consequence of the Gaussian 
estimates in Proposition 4.3 it turns out that these estimates carry over to any 
L P (L P ) space, where L p := L p (pL a ), provided that p G (1, oo). We will thus call I, k, 
and fi in p e d^dfw Calderon-Zygmund exponents whenever 

|/3|)g {(0,1,0), (0,0,1), (1,0,2)}. 

For these exponents, the kernels p £ dfdfG satisfy so-called Calderon-Zygmund can¬ 
cellation conditions, which in turn yield the L P (L P ) estimates. Going one step 
further, Muckenhoupt theory allows for dropping the weights if p < oo is chosen suf¬ 
ficiently large. We thus deduce maximal regularity estimates in regular L p = L P (L P ) 
spaces. 

We refrain from establishing the cancellation conditions in this paper, which would 
be a straightforward but tedious calculation based on the Gaussian estimates of 
Proposition 4.3, We refer the interested reader to Kienzler’s work [34], specifically 


Corollary 3.18, therein. 


Lemma 4.13. Let w be the solution to the initial value problem 4.33, (4.34). Then 
for any p > max {l, yyy} it holds 

\\d t w\\ L p + ||V?u||lp + ||pV 2 H| LP 

and 


lp 


\\d t Vw\\ LP + ||V 2 w;||lp + Wp^Mlp < ||V/|| LP . 
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Proof. We first notice that Lemma [4.6| ensures that the mappings 

L*(Ll) 3 / pftfdZw e r(Ll) 


are continuous if £, k, and (3 are Calderon-Zygmund exponents. The correspond¬ 
ing kernels, p e d^d^G, satisfy the Calderon-Zygmund cancellation condition, cf. [33, 
Corollary 3.18], and thus, maximal regularity theory applies and yields 

||A‘3fHI mm < ll/lk-M). 


Furthermore, the Muckenhoupt theory allows for dropping the weights. More pre¬ 
cisely, p e ~ a is a p- Muckenhoupt weight with respect to p a if and only if — 1 < 9 < 
p(a + 1 ) — 1 . We can thus translate the maximal regularity theory from L P (L P ) to 
L p (L p (p 9 ~ a pa)) = L p (Lq). In particular, for 9 — 0, it holds 

< II/IIl- (4.35) 


provided that p > max {l, 

The second estimate is obtained from (4.35) by differentiating the equation and 


iteration. Thanks to strict parabolicity, this amounts to a standard exercise in the 
interior of the unit ball. Close to the boundary, however, we have to argue more 
carefully. In fact, after a suitable localization, the techniques developed in the proof 
of Lemma [4.8| are applicable. We claim that 


H/afafvHiw < II/IMjkb) + l|v/|U„ 


and leave the details to the reader. Let now h(t) := |-Bi(0)| -1 / f(t,z)dz and 
H(t ) := f* h(r ) dr. Obviously, w — H is a solution to the equation with inho¬ 
mogeneity / — h, and the previous estimate applies. Invoking Poincare’s inequality 
11 / — h\\LP < ||V/||lp and observing that d — H) = d^d^Ww for any choice 
of k and j3 yields the second statement of the lemma. ■ 

In the proof of Proposition |4.4| it will be convenient to decompose / = ?// + (1 — p)f 
where p denotes a cut-off function which is constantly 1 in the cylinder Q^.(zq) and 
vanishes outside the larger cylinder 

Qr(z 0) := X 5 2r(-0)- 


First, we provide a “diagonal” L p estimate. 


Lemma 4.14. Let w be a solution to the initial value problem (4.33), (4.34). Sup¬ 
pose that spt / C Q^{z o) for some z 0 G -£>i(0) and 0 < r < V27 Then for all 
Calderon-Zygmund exponents l, k and /3 and p > max {l, it holds 


r 2 \Q d r (zo)\ p\\p e d^d^w\\ LP{ Qd {zo)) < ||/||y ( p). 
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Proof. We deduce from Lemma [4.13| the estimate 

llA*afv w ||„ (0 <, w) < ||v/||„ (ft(w)) . 

Let {Qf.(zi)}i e i be a finite cover of Q d (z o) such that r; ~ r and \Q d .{zf)\ < 
|Qr(^o)l- Then 


liv/ll 


LP(Qd(z 0 )) ^ 


< 5L ii v /iu-(os ! 4 

i£l 


(*i)) 


< i- 

rsj 2 




ie/ 


< ^i^(^)hi/iiF (P ). 

Notice that |®(« 0 )| = |r J |.B| r (^o) | < r 2 |^(^ 0 )| = |<5r(-o)| by Lemma 
proves the lemma. 


4.3 


This 


In the case (£, k,\/3\) = (0, 0,1), if \Jp{zf) r, the estimate of Lemma 4.14 is not 
strong enough to obtain a Carleson measure bound. The stronger estimate will be 
derived in the following lemma. Simultaneously, we establish a pointwise estimate 
on Vic. 


Lemma 4.15. Let w be a solution to the initial value problem (4.33), (4.34). 
Suppose that spt / C Qf(zo) for some zq G -£>i(0) and 0 < r < y/2. Let p > 
max {iV + 1, ~p[ } ■ Then, for any t G [0, r 2 ] it holds 

| Vtc(f, Zq)\ < \\f\\ Y (p)- 

Moreover, if ^p(zo) 3> r and p > max {N + 2, it holds 

r(r+ y/p{z o)) \Q$(zo)\~* || SJ 2 w\\ LP{Q a {zo)) < ||/||y (p ). 

Proof. It is convenient to deduce the two statements from the following estimates 
on w and Vic, respectively: Let / G L p . Suppose that w is a solution to (4.33), 
(4.34) with / replaced by /. If p > N + 1, then 

Ht,z)\<r 2 \Q d r (z)rHf\\ LP (4.36) 


for all (t, z) G [0, r 2 ] x L?i(0), and if p > N + 2 and \/ p{z) > r then 

ry/p(zj\Vw(t,z)\ <r 2 \Q d r {z)\~p\ 


(4.37) 


for all (t,z) G [0,r 2 ] x Lu(0). 

Indeed, differentiation of ( |4.33 ) with respect to the i-th coordinate gives the equation 
dfdiW — p _cr V ■ (p CT+1 V9jU;) = dif + ZiAw + (a + 1) (d,w + z ■ Vditv). 
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Let / be the right-hand side of this identity. Applying (4.36) to w — diW and 
su mming over i yields 

|Vio(M)l < r 2 \Ql(z)\-i (||V/||i, + ||Vtu||„ + ||V 2 <o|M, 

using that / is supported 


for all (t, z ) G [0, r 2 ] x i?i(0). With the help of Lemma 
in Q d (z o) and choosing z = Zq, the latter turns into 


4.13 


|Vw(t,^ 0 )| <r 2 \Q d (z 0 )\ * 


W^P(Qd(z 0 ))- 


Arguing as in the proof of the previous lemma, the right-hand side is controlled by 
11./' 11 Y(p) i which shows the first statement of the lemma. 

In a very similar way, we deduce from ( |4.37[ ) the estimate 

r^fp{z)\V 2 w(t,z)\ < \Q$(z)\-p\Q$(zo)\>\\f\\ Y (p), 


for all (t,z) G [0,r 2 ] x B d (z 0 ). In view of Lemma 4.2, it is p(z) ~ p(z 0 ) for all 
z G B d {z 0 ) provided that p(z 0 ) r. Hence \Q d (z)\ ~ \Q d (z 0 )\ and integrating 
the previous estimate over Q d (zo) yields the second statement of the lemma. 

It remains thus to prove (4.36) and (4.37). Both estimates can be established simul¬ 
taneously. We use the heat kernel to represent the solution, cf. Lemma 4.12 , and 
estimate with Holder’s inequality 

\d^w(t jZ )\ < || d%G(t, - ,z, 


for all (t,z) G [0,r 2 ] x L>i(0), and where ^ ^ = 1 as usual and /3 G such 

that \j3\ G {0,1}. The estimates (4.36) and (4.37) follow thus from appropriate L q 
estimates on the heat kernel. First, recall the Gaussian estimate from Proposition 


\dPG{t,t',z,z')\ < Vr 1/31 (Vr+^/p(z)) ' \B d ^{z)\ a x e 


?d f „\\-l -Cd(z,z') 2 /r 


(4.38) 


where we have set r = t — t! and we have used Lemma 4.4 to replace the ball centered 


at z’ by the ball centered at 0 . We let | B d ^(z)j be a cover finite cover of 

-Bx(O) and estimate 




Invoking Lemma 
for some k > 0. 


4.3 


3 = 1 

< yy«c«-D s | B^z) 

3 = 1 

j K \B d 


e-gC^z^/r 


qa- 


ieiice 


we compute \B d jV ^{z)\ qc < j K \B^(z)\ qa ~ j K \B d ^(z)\ l q \B d ^(z)\ 


e - qCdM2 lr d p q(J {z') < ^ j K e~ gC ( j ~ 1)2 \B^(z)\ 1 - q \B , ^ r (z)\ q a . 

je N 
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Because the series is convergent, an integration of (4.38) over (0,r 2 ) x -Bi(O) yields 


I J \S^G(t,t',z,z , )\ q diJ. w (z')dt' < J 


T~ m I VT 


_\ —q\(3\ 

+ VM) iBywi 1 -’*. 


Once again, we distinguish the cases where z is relatively close to the boundary 
and the case where it is not. In the first case, if a/ p(z) < r, we choose (3 = 0 and 


estimate with the help of Lemma 4.1 


/ | B^z)^ dr < / r^ N dr <r 2 ^ N+2 <r 2 \B?(z)\ 1 ~ q , 

Jo Jo 

provided that p > N + 1. In the second case, if \Jp(z) > r, we have 

, -m 


< 


Vr qm (y/r+ y/p(zfj 9 | B d ^(z)\ l q dr 


V 7 /^) 


- q \/3\+(l-q)N 




< 

r^j 


^-^- q \p\+^-q) N r _ q \ p \ Hl - q)N+2 

provided that (2 — \/3\)p > N + 2. Therefore, in either case, if p > N + 1 it holds 
11 ^ 0 , -,z, Oil L 9 (( 0 ,r2);L^) ^ rlq \ B *(z)fi- 1 ~ r 2 \Q d r (z)\~p, 


which implies (4.36). The estimate (4.37) follows analogously with |/3| = 1 and 
p > N + 2. ■ 

We come now to the “off-diagonal” estimates. We define 

^ „2 


Qr(z o) := ( yr ,r j x Bi(z 0 ), 


so that Qr(zo) c Qr(zo) C <w0o)- 


Lemma 4.16. Let u> &e the solution to the initial value problem (4.33), (4.34). 
Suppose that spt / cj0,r 2 ] x L?i(0) \ Qr{z 0 ) for some z 0 G £?i(0) and 0 < r < a/2. 
Then for all (f, z) G Qr(z 0 ) an d P > max {1, yyr} ^ holds 

|Vie(f, z)| + r 2 |(9 t Vw(t, z)l + r + \fp(z) j | V 2 re(t, z)| + r 2 p(^)| V 3 re(t, z)| 

~ I|/||v(p)- 

Proof. We start the proof with an elementary auxiliary estimate whose proof we 
postpone until later: We claim that there exists a constant C > 0 such that 


yjl _ p ° e -Cd{z,z')/(t-t') < r -e e -Cd(z,z')/r 


(4.39) 
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for all (t,z) G Qr(z o) and (t', E) G [0,r 2 ] x L>i(0) \ Q^(zo). In the following C will 
denote a nniversal constant whose value may change from line to line. 

For all (t, z) G Q^{ z o) if is /(f, z) — 0 and thus 


KSfw(M) | 


y/r 2k ^ ( v r + v I j e cd ^' z ')‘ t /r 


< 

nsj 


\f(t - r, z')| d/v(z')dr 


by Duhamel’s principle (Lemma 4.12) and the Gaussian estimates (Proposition |4.3[ ) . 
The estimate ( |4.39| ) and the monotonicity of the function s H* ^ show that we 


can replace yfr by r in the above estimate. Thanks to Lemma 4.4 


we can moreover 


substitute balls centered at z' by balls centered at 0 and vice versa. We thus have 


r 2fc+|/3|-l . r 


< 


+ 

r e -Cd(z,z')/r 
r\ B r(z)\a 


1-1 


\d*dPw(t,z)\ 


r + a /p(z')) \f(t',z')\dp, a (z')dt'. 


Let {Br(zi)} 1<i<n be a family of balls covering B i(0). Observing that d(z,Zi ) < 
d(z, z') +r for any z! G B^(zi), we further estimate the expression on the right-hand 
side by 


£• 

, t=1 


—Cd(z,zi)/r 


e ~Cd(z,z)/r rr 2 r- 

SS r\B*(z)\ a J 0 j B d (i) 


r + y/p ) 1 I/I dp a (z')dt'. 


The sum is bounded uniformly in r and z. Moreover, thanks to Lemma 4.4 


we can 


replace z by z in the multiplicative factor in front of the integral. Hence, to prove 
the statement of the lemma, we have to show that 


\B d (z)\ + •/'llL 1 ((0,r 2 )xBd(5)) ~ ll/llx(p)- (4.40) 


For that purpose, we cover the domain of integration with a countable number of 
intrinsic cylinders. More precisely, with Tj = y/3/2 "V for j G N 0 , we hnd a hnite 

family of balls Bf.(zij) \ with Zij G B^{z) such that 

l 3 J 1 <j<rij 


n-j 


(4.41) 


uniformly in j and r. (The latter follows from a non-Euclidean version of Vitali’s 
covering lemma, cf. [36, Lemma 2.2.2]). Then (0 , r 2 ) x B^(z) C U,- e N 0 Ui<i< n3 .(3r. (z^) 
and thus, via Holder’s inequality 

_ n i 

l|p' T (r + v / P)“ 1 /ll Ll((0 , r 2 )xB^(z)) ~ E E n x G + v ^) _1 II Li{Qd, (.zij )) \\f\\LP(Qd, ( Zij )), 

JGN 0 i=l 
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where - + - = 1. 
p q 

Bi 


It follows from Lemma 


4.2 


that ( rj + y/p) 1 < {rj + ^p(zij)) 


-l 


m 




Therefore, 


II P a ( r j + Vp) II LHQiUij)) 


< 


r i+\/p(%) I QrM 


-1 


Hj)\qv 


\fp( Z ij) 




Notice that the measure p qa is finite only if qa > —1, or equivalently, p > ^rry. 


Combining the previous two estimates and using (4.41) and the convergence of the 


geometric series yields (4.40). 


We finally turn to the proof of (4.39). we first consider the case where z' G B d (zo). 
We then must have t' / (^, r 2 ), and thus r = t — t' > /r 2 . Then 

e ~cd{z,z ') 2 /t < y/r~ e < r~ e < r ~ e e ~ Cd ^ z,z '^ r , 

because d(z,z') < d(z,zo) +d(zo,z') < §r. In the complementary case, where 
z' / B d (zo). it holds that \ < d(z, z') because z E B d (z 0 |^ Now, a small computation 

shows that the function r i—>■ e~ Cd ^ z,z '^ 2 / T y/r 6 is increasing as long as 0 < r < 
d(z, z') 2 . Since r < r 2 < d(z, z') 2 , it thus holds that 

^~ e e ~Cd{z,z') 2 /T < r ~e e ~Cd(z,z') 2 /V 2 <- r ~e e ~Cd{z,z')/r 


This proves (4.39). 


We are finally in the position to derive the Caleson-measure estimate for the initial 
value problem (|4.33|), (14.341). 


Proposition 4.4. Let w be a solution to the initial value problem (4.33), (4.34). 
Suppose that p > max {iV + 2, j. Then 


MU(p) + ll w llLi P + ^ ll/l|y( P )- 


Proof. We start with the estimates for small times, i.e., we assume that r < a/2. 
There is no loss of generality to assume that f(t , •) = 0 for t > r 2 . We fix (t, £o) G 
(0,r 2 ) x B\ (0) let rj denote a cut-off function such that rj — 1 in Qf(zo) and rj = 0 
outside Q d (zo). We decompose f — f\ + f '2 with /1 = rjf and / 2 = (1 —17)/, so that 
spt /1 C Q d {zo ) and spt / 2 C [0, r 2 ] x B \Q d (zo). Let w\ and w 2 denote the solutions 
to the inhomogeneities /1 and / 2 , respectively. By the linearity of the equation, it 

4 Here we assume that the triangle inequality holds for d for notational convenience. In fact, 
to guarantee a lower bound d(z, z') > r in this case, we would have to consider intrinsic balls of 
smaller radius in the definition of Qf(zo). 
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holds w = Wi + w 2 . We deduce the control of w\ from Lemmas 4.14 and 4.15 and 


of w 2 from Lemma 4.16 For instance, for all t G [0,r 2 ], 

\Vw(t,z 0 )\ < \Vwi(t, Zo)\ + \Vw 2 (t,z 0 )\ < ||/i||y(p) + \\f 2 \\Y(p)- 
The latter is controlled by ||/||y( p ) because 

|V?7| < r ” 1 (r + v / p5o)) • 


We thus have \\7w(t,z 0 )\ < ||/||y( P ) for all (t,z 0 ) G [0,r 2 ] x L?i(0). The right-hand 
side being independent of r and zo, this proves the desired Lipschitz bound for w 
for small times t G [0,2], The Carleson measure estimates in this time interval are 
derived similarly. 

ft remains to derive the estimates for large times, T > 1. By superposition, in 
view of the above estimates and Proposition |4.2 we can thus assume that / = 0 in 
[0,1) x l?i(0). We let x — X(t,t+ i) and decompose f — fi + f 2 with fi = xf and 
/ 2 = (1 — x)f- L e f ^1 and w 2 denote the solutions to the inhomogeneities fi and 
/ 2 , respectively. 

From Lemma [4. 13 [ we deduce that 

||^tVu'i||z / p(Q( T )) + ||V 2 u;i||lp(q(t)) + ||pV 3 -u;i||lp(q(t)) < ||V/i||lp(q(t))- (4.42) 

To gain control on the Lipschitz norm of w 1 , we invoke the Sobolev embedding 
theorem, which, if p > N reads 


IICIIl°°(q(t)) ^ IICIIlp(q(t)) + ||^CI|lp(q(t)) + II veil lp(q(t))- 


Choosing ( = Vu)\ and using a Poincare estimate in time with Wwi(T) 
yields 

||Vu;i||loo(q( T )) < ||V/i|| L p(q(t)) < ||/i||y( P )- 


0, (4.42) 


In order to derive the analogous estimates for w 2 , we use the exponential decay of 
the Gaussian in Proposition 4.3 to estimate 







\f 2 {t',z')\ dn a (z')dt', 


for all (t, z) G (T + |, T) x i?i(0). We employ Holder’s inequality and use that ^i qa 
is a hnite measure if its dual p > —C-. Then 

1 cr+1 


1 T 1 / rn +1 \ 1/9 

\S?d?w 2 (t, 2)1 < v / dt ' ||/ 2 |m (Q(fc)) . 

n =1 \ Jn ' 

The sum is estimated by e Xlt , and we obtain thus control in the Lipschitz norm 
and the L p norm in the cylinder (T + |,T) x fi^O). Arguing as before, we obtain 
control of w = W\ + w 2 in that cylinder because ||/i||y( p ) < ||/||y( P ), and the bound 
is independent of the decomposition. The desired statement easily follows. ■ 
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5 The nonlinear problem 


In this section, we turn to the proof of Theorem 2.1, that is, we will show well 


posedness of the nonlinear perturbation equation (]2.5|) in a neighborhood of the 
constant ta* 


m 


= 0 and establish analytic dependence of solutions on the initial data 
that neighborhood. For convenience, we rewrite ( |2.5[ ) as: 

d t w - p~ a V ■ (,o u+l Ww) = /3pF(w, Xiu) - p~ a V ■ (,o a+1 zF(w , Vw)) , (5.1) 


where f3 = IV + 2a + 1 and 


F(q,p) = 


IpP 


1 + q + z ■ p 


for (g,p) G R x R A . For abbreviation, we will sometimes denote the right-hand side 


of (5.1) by f(w). Well-posedness for this equation will be derived by combining the 
linear theory of the previous section with a fixed point argument, that relies on the 
observation that the nonlinearity is essentially quadratic in a neighborhood of the 
constant re* = 0. However, a direct application of the linear theory presents us with a 
technical problem: The (semi-)norm on the solution space X(p) in Section [4] contains 
only the homogeneous part of the Lipschitz norm, whereas the nonlinearity F above 
demands control of the full C 0,1 norm to prevent the denominator from degenerating. 
We will face a second problem later in Section [6j Theorem |2.1 shows that the 
nonlinear problem (5.1) generates a C 0,1 semi-flow locally in a neighborhood of the 
constant w* = 0. However, in order to fit into the dynamical systems framework 
needed to establish invariant manifold theorems, it is necessary to extend the local 
flow into a global one. 

We will tackle both problems simultaneously by considering a truncated version of 


(5.1) first. More precisely, we introduce a smooth function that cuts the nonlinear 
terms off at points where the solution or its gradient is too large. The resulting 
equation is linear at such points. If, however, a solution is globally sufficiently 
small in C 0,1 , then the cut-off is inactive and the truncated equation coincides with 
the original one. The latter is guaranteed by a smallness condition on the initial 


datum thanks to the a priori estimates from Theorem 4.1 and a (weak) comparison 


principle, cf. Lemmas 5T and 522 below. 

With regard to the invariant manifold theorem that will be derived in Section [6j 
working with the truncated equation instead of (5.1) has one more advantage: The 
truncated equation generates not only a semi-flow in C 0,1 , but also in L We prefer 
to exploit the latter as it closely connects to the spectral analysis conducted in B3 
Moreover, many of the concepts used in the derivation of the theorem are easier 
explained and formulated in the Hilbert space setting. We will eventually transfer 
the results from L 2 to C 0,1 with the help of smoothing estimates. 


This section is organized as follows: In Subsection 5.1, we introduce and study the 
truncated problem. We prove well-posedness in L 2 and establish a priori estimates 


in smoother spaces. Subsection 5.2 contains the proof of well-posedness of the full 


nonlinear perturbation equation, Theorem 2.1 
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5.1 The truncated problem 

We let f) : [0, oo) —> [0,1] be a smooth cut-off function that is supported on [0, 2) 
and constantly one on [0,1]. Let 0 < e, S < 1 be fixed such that a/2 (e + 5) < 1. We 
then set for (q,p) G R x R v 


/ \ M r w / \ p I 2 

VeM’P) = T i\^) r i\—r 


and with that F £ j = Ve,sF. The truncated equation is 

d t w - p- ff V ■ (p ff+1 Vw) = (3pF £) s(w, Vw) - p- ff V ■ (p a+1 zF £ ,s{w, Vw)) . (5.2) 

For abbreviation, we denote the right-hand side by f £t s(w). It is clear that this 


equation is equivalent to (5.1) as long as \w\ < 6 and |Vw| < £. 


Notice that in (5.2), the cut-off rj £j s acts on both w and Via, so that the nonlinearity 
F Ej s(w, Vw) is globally defined. Conveniently, the cut-off not only allows for a global 


c o,i 

semi-flow as a consequence of the estimates derived in Theorem |4.l[ It also has 
a regularizing effect as now a Hilbert space theory becomes applicable. Via a fixed 
point argument, it yields the existence of a global L 2 semi-flow. It seems hence 
natural to first construct invariant manifolds for the flow for the truncated equation 
in the Hilbert space setting and then carry the results over to the original equation 
for the C 0,1 semi-flow. This will be done in Section H] below. 

We finally remark that for the derivation of Theorem |2.1[ the study of the truncated 


equation is redundant and Theorem 2.1 can be established directly using Carleson 
measure estimates and a fixed point argument exploiting the quadratic behavior of 
F(q,p) in p near the origin. For this, however, in addition to the Carleson measures 
in Theorem |4.1[ which are formulated for the gradient of w, we would have to 


establish the analogous bounds for w itself to gain control over the full C 0,1 norm 
in terms of the initial datum. In this regard, the situation considered in this paper 
is slightly more complex than the one studied in [34] . 

We now state and prove 

Theorem 5.1. There exists e 0 > 0 such that for every 0 < e < e 0 and 0 < 5 < 1 
with v^2 (e + 5) < 1 the following holds: For every g G L 2 , there exists a unique 
solution w to the truncated equation (5.2) with initial datum g. Moreover, w G 

L“((0,oo)|Lj)nL 2 ((0,«D);//i +1 ). 


i 2 (iGi). g G L 2 , and T > 0, the initial value problem for 


Proof. For every F G 
the linear equation 


d t w - p“ CT V • (p CT+1 Vw) = (3pF - p- ff V ■ p a+1 zF J (5.3) 

has a unique solution in the time interval (0, T ). It satisfies the a priori estimate 


\W L‘ 


= (L2) + ||Vw|| l2 ( L 2 +i ) < C T (||T|| L 2 (L 2 +i) + ||p||z,2 j , (5.4) 
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where Ct is a generic constant that depends on time T. (Here and in the follow¬ 
ing discussion, the norm dependences on the time interval (0, T) are suppressed.) 
Existence follows from standard Hilbert space methods and ( |5.4 ) can be derived by 
testing the equation with p a w — we omit the details. 


To prove well-posedness for the nonlinear equation (5.2) for any initial datum in L 
we invoke a fixed point argument. Given a function w G L°°(L' 2 ) fl L 2 (H^ +1 ), we 


2 

(7 ? 


denote by w = w(w,g) the solution to ( 5.3| ) with F = F £) s(w,X7iu) and with initial 
value g. Because \F £jS (q,p)\ < e\p\, we have ||^||i, 2 (i 2 } < e||Vu;|| £ 2 (L 2 } , and thus 


by (5.4) 


+ C T 


Ll 


IHU°°(n|) + I|Vw||i, 2 (L 2 +i) < G t£ ||Vu;|| L 2 (£ 2 +i ) 

for some universal (but time-dependent) constant Ct > 0. Moreover, because 
\F E ,s{qi,Pi) -F £j5 (q 2 ,p 2 )\ < £ (|<?i ~ ^1 + \Pi - P 2 I) — as the reader may straightfor¬ 
wardly compute — applying (5.4) again, we obtain for two solutions w\ = w(wi,g) 
and W 2 = w(w 2 ,g) that 

\\w 2 - wi\\l°°(l%) + II Vu ) 2 - Vwi|| L 2 (L 2 +i) 

< C T e (\\w 2 - wi||l°°(l 2) + ||Vw 2 - VH| L 2 (L 2 +i) ) • 

Hence, for e sufficiently small, the solution map w 1 —y w(w,g) is a contraction in 
L-(Ll) n L 2 (Hl +1 ), and thus, for such e, there exists a unique local solution in 
L°°(Ll) fl L 2 (Hl +l ) to (5.2) with initial datum g G L 2 . 


To show the existence of a global solution, it suffices to derive a bound on solutions 


of (5.2) that is uniform in T. For this purpose, we test the equation with w — c 
where c > \/2 5 is a fixed constant. Then 


d 1 f. 

-- n w -c) + 


\Vw\ 2 dfJLv+i 


= (3 j (w — c)F £j s(w, Vw) dp a+ i + J (z ■ Vtv)F £} s(w, Viu) dp a+ 1 . 

Since 0 < F £j §(w, Vw) < e|Vtu|, the second term on the right can be easily absorbed 
into the left-hand side provided that £ is chosen sufficiently small. Moreover, because 
\w\ < \/25 in the support of F £t $(w,X7w), the first term on the right-hand side is 
non-positive: 

/ (w - c)F £t5 (w , Vw) dpa+i <(S — c) / F £ j(w, Vw) dp a+1 < 0. 


Thus 


^ [(w - c) 2 dp a + f | Vrc| 2 dp a+ 1 < 0, 


for some universal constant C > 0. Integration yields the desired uniform bound. 
This concludes the proof of Theorem 5.1| 


Our Erst higher regularity estimate is for Lipschitz initial data. 
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Lemma 5.1. There exists £o > 0 such that for every 0 < £ < £o and 0 < 5 < 1 with 
y/2 (£ + 5) < 1 the following holds: If w is the solution to (5.2) with initial datum 
g € C 0,1 , then w is smooth and, for any p £ [ 1 , oo) it holds 


IMIa-(p) + IMIlip ~ Iblkip- 

Proof. It is enough to prove the estimate for p > max jiV + 2, W_}. The extension 
to the full range p G [1, oo) then follows by Holder’s inequality. If £ 0 , £ i an d 6 are 
as in the assumption of Theorem |4.1[ then we have the estimate 

IM|x(p) + \\ W \\u P < \\fe,s(w)\\ Y {p) + 111?11Lip- (5.5) 

Therefore, the statement of the lemma follows provided that 

||/ e ,«HI|r(p) < ^ (IMU(p) + ||w|| L ip) , (5.6) 

and upon choosing £ 0 (and thus e) smaller if necessary. 

We start noticing that \f e ,s\ ^ p|VF| + p|F||Vw| +£~ 1 p\F\\'V 2 w\ + |F| and |V/ £ ^| < 
|pV 2 F| + £ _1 p|V.F| + |VF| +£~ 2 p\F\ + £ _1 |F| in the support of the cut-off function. 
(The computations leading to these and the following estimates are straightforward 
but tedious. We omit the details.) Hence, by chain rule and the definition of F, 

\fe,s\ ;$ £p\^ 2 w\ + £|Vw| 


and 

|V fe t s\ < e\ Vw| + e\ V 2 w| + p|V 2 u >| 2 + £p|V 3 w| 

in the support of r} £t s- In the estimates of f e j and Vf Si s, the terms that deserve 
a special consideration are those which have second order derivatives, i.e., £p|V 2 w| 
and e\ V 2 w| 2 . We restrict our attention to the estimates for the time interval (0,1). 
First, an application of Lemma |4.2| yields 

-- 7 ===||pV 2 ry|| L p ( Qd (2)) <r(r+ y/p[z)) || S7 2 w\\ LP{Qd{z)) . 

r + Vp( z ) V ' 

For the quadratic term, we use the interpolation inequality || V£|| b p ~ IICIU°°IIv 2 ciUg, 

L/p P 

cf. [33, Proposition 2.18], and obtain 


|V W | || LP(Qd( 2 )n S pt r/ e s) ~ IIVHI L°°(spt r/ £ s) ||pV w|| lp (Qd( 2 )n S pt r) Si s)' 


Then (5.6) follows because ||Vw||l°c < £ in the support of r) £> , 5 . 
Since w is Lipschitz, (5.2) can be seen as the parabolic equation 

d t w - p- CT V • (p CT+1 ( 1 - H)Vw) = f3pF E , s (w, Wiu ), 


where A given by 


A = g £ ,s (w,Vw): 


Vw 


1 + w + z ■ Vw 

is bounded and 1 — A is elliptic. Hence, higher regularity can be obtained from 
regularity theory for parabolic equations with rough coefficients, see, for instance, 
Theorem 5.6.1] for a theory in the degenerate setting. ■ 


56 










We finally establish the following weak comparison principle for the nonlinear equa¬ 
tion (5.2). 


Lemma 5.2. Let £q, e and S be as in Lemma 5.1. Let w be the solution to (5.2) 
with initial datum g G C 0,1 . Suppose that a < g(z ) < b for some a, & G R and all 
z G -Bi(O). Then 


a < w(t, z) < b for all (t, z) G (0, oo) x £>i(0). 


In particular, 


I W || J^OO < f 


Proof. By Lemma 5.1, w is smooth, so that we can give a pointwise argument. 
We only prove the lower bound, the upper bound is established analogously. Let 
7 > 0 be arbitrarily given. We will show that a — 7 < w(t,z ) for all t > 0 and 
z G Ri(0). Since 7 > 0 was arbitrary, the result follows. We argue by contradiction. 
Suppose that there exists i* > 0 and z* G -£>i(0) such that a — 7 = w(t*,z*) and 
a — 7 < w(t, z ) for all t < i* and z G -£>i(0). In the case where z * lies in the interior 
of the domain, z * G L>i(0), then d t w(t*,z *) < 0, Vw(t*,z*) = 0 and X7 2 w(t*, z*) 
is positive definite. It immediately follows that ( d t w — p _<7 V • (p ff+ 1 Vty)) < 

0 while fe,s( w ,)(t*, z *) =0 — a contradiction. Otherwise, if z* G <9Ri(0), then 
d t w(t*,z *) < 0 and z * ■ Vrc(t*,£*) < 0. We deduce that 

(d t w - p~ a \7 ■ (p ff+ 1 Vw)) (£*, z*) = d t w(t *, 2 :*) + (a + 1)^* • 2 *) < 0, 

which contradicts f(w)(t *, 2 :*) = (cr + l)F(ta, Ww)(t *, z*) >0. ■ 


5.2 The full problem. Proof of Theorem 2.1 


It is convenient to make the following definition: For e, 5 > 0, we consider 
B e ,8 '■= {w G C 0 ’ 1 : ||w||i/x= < 5, 11re11Lip < s} ■ 

We are finally turn to the 


Proof of Theorem 2.1. Existence of a Lipschitz solution for sufficiently small initial 


data follows immediately from the theory for the truncated equation (5.2): Theorem 


5.1| provides us with a unique solution for the truncated equation, which is also a 


solution to the original equation thanks to Lemmas 5.1 and 5.2 It is clear that this 


solution is unique by the uniqueness for (5.2). 


We shall show that solutions depend analytically on the initial datum. For this 
purpose, we restrict our attention to finite time intervals (0, T), but suppress the 
dependence on T in the following. We first notice that both 


and 


B e ,5 F(w,Vw) G L 2 (L 2 +1 ), 

X(p) n B £)& 3 f(w ,) G Yip) 
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we deduce from Lemma 4.5 and Theorem 4.1 


are analytic mappings. Moreover, because of the elementary estimates 

IlClUl + IlCIILip and IICIUj < IKII 

that 

< W f \\l^lI +1 ) + \\f\W(p) + 


< 


|ic||co,i 


c°a- 


Therefore, if w — w(w,g ) denotes the solution to the linear perturbation equation 
with / = f(w) for some w £ X(p) D B £y § and initial datum g £ B £ j, then w : 
X(p ) fl B e j x B e ^ -£ X(p) D C 0,1 is analytic. Hence, also S : X(p) fl B e< s x B e> s —> 
X(p) fl C 0,1 with E(w, g) = w — w(w, g ) is analytic. It moreover satisfies S(0, 0) = 0 
and £^5(0,0) = 1 , and thus, the analytic implicit function theorem shows (upon 
decreasing S and e if necessary) that the unique solution w = w(g) of E(w,g) = 0 
depends analytically on g. However, by the definition of 5, w is the unique solution 
to the nonlinear perturbation equation (5.1) with initial datum g. 

The proof (2.6) proceeds along the same lines as Kienzler’s proof of the analogous 
statement in [33, Theorem 4.1]. Only instead of distinguishing between “horizontal” 
and “vertical” directions, we should argue as in the proof of Lemma T 8 , that is, we 
rather distinguish between “angular” and “radial” directions. ■ 


58 







6 Invariant manifolds 


In this section, we turn to the proof of Theorem £2 

The operator obtained from linearizing the truncated equation in the Hilbert space 
L 2 around the constant function w* = 0 differs from the linear operator studied 
in m only in the underlying Hilbert space (and a multiplicative factor, that we 
neglect in the following discussion): As in [47] . the action of the linear operator £ 
on a smooth function w is given by 


\7p ■ S/wdp a+ i = / LpCw dp a for all cp £ C oo (H 1 ( 0 )), 


( 6 . 1 ) 


i.e., Cw = —pAw + (cr + 1 )z ■ Vw for smooth w. The role of boundary conditions 
for this elliptic problem is discussed in Remarks 4.1 and 4.2 in Section [4] above. 

In the present paper, we choose the Lebesgue space L 2 as the underlying Hilbert 
space, while in m, the author diagonalized the operator in the Sobolev space 
,_i/R. To distinguish these two operators, we will in the following denote by Tt 
the operator on H^ +1 /R. 

We shall be more specific with regard to the underlying functional analytical con¬ 


cepts. The operator defined via (6.1) on C' oo (I? 1 (0)) is a densely defined 


nonnegative symmetric operator both on L 2 and on Hf +1 /R (cf. [47, Appendix]). 
Hence, £|c°°(spo)) * s c l° sa ble in both Hilbert spaces and we denote its closure in L 2 
by £, and its closure in Hf +1 /R by TL. Since both operators are defined through 
the closure of a quadratic form, £ and TL are Friedrichs self-adjoint extensions (cf. 
[ 461 Chapter X.3]) of £|c , oo(^y) in L 2 a and H\ + i/R, respectively. 

Comparing functions in HA ±/R and L 2 reveals some information on the spectrum. 
In either case, eigenfunction are smooth and lie in any Hilbert space (T+k . Hence, 
even though £ is set in a larger space, we expect that 0 is the only new eigenvalue 
in the discrete spectrum, i.e., 


£(£) = E{n) U{0}. 

The emergence of the new zero eigenvalue stems from the fact that L 2 contains the 
constant (eigen-)functions, which were rnodded out in the quotient space H^ +1 /R. 
This observation is the content of the proposition that follows. The computational 
results are taken from m and are appropriately rescaled. 

Proposition 6.1 (The spectrum of the linearized operator). The spectrum of £ 
consists entirely of isolated eigenvalues with finite multiplicity. They are given by 

\ £k = ( a + 1)(£ + 2k) + k{2k + 2t + N - 2), 

where (£, k ) £ N 0 xN 0 if N > 2 and (£, k ) e{0,l}xN o */7 = 1. The corresponding 
eigenfunctions are polynomials of the form 


ipoio(z) 
‘fink (t) 


1 , 

F ^—k, a + £ + — + k;£ + —; |^| 2 ^ Yg n 
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where n E {1,..., Ng} with Ng — 1 if £ = 0 or £ = N = 1 and Ng = (jV+£ £ ,^ ( ^^, 2f 2 - > 
else. Here, F(a,b;c;z ) is a hypergeometric function and in the case N > 2, Yg n is 
a spherical harmonic corresponding to the eigenvalue £(£ + N — 2) of A s n-i with 
multiplicity Ng. Otherwise, if N = 1, it is Ya(±l) = (±1) £ . 

The literature on hypergeometric functions and spherical harmonics is vast, see, e.g., 

mmmm- 

Hypergeometric functions F(a, b ; c; z) can be defined as the power series 


F(a, b ; c; a;) 


i , v" ( a )j( 6 )j ^- 


for any a, b, c, x E R and c is not a non-positive integer. The definition involves 
the extended factorials or Pochhammer symbols 


(s)j — s(s + 1) ... (s + j — 1), for j G N, s E R. 


The power series is convergent if |x| < 1. In the situation a hand, hypergeometric 
functions come into play in the study of the radial part of the eigenvalue equation. 
The latter is a second order Fuchsian ODE with three regular singular points that 
can be transformed into the hypergeometric differential equation. In [57], the author 
solved this equation and selected those solutions as eigenfunctions that lie in the 
domain of the linear operator. 

Plugging the value a = —k into the definition of the hypergeometric functions, 


we see that in Proposition 6T the hypergeometric functions reduce to polynomials 
of degree k in x = \z\ 2 . Recalling that spherical harmonics Yg n are homogeneous 
harmonic polynomials of degree £ (defined on the unit sphere), we thus observe that 
the eigenfunction ifg n k is a polynomial of degree £ + 2k. 

At this point, we shall briefly compare our spectrum with that found by Denzler and 
McCann in EH for the linearized fast diffusion equation. Since the solutions to the 
fast diffusion equation are positive instantaneously, the linearized operator is defined 
on all of R iV , and thus, due to spatially algebraic decay of solutions, the operator 
has non-compact resolvents. It follows that the spectrum is non-discrete. To be 
more specific, Denzler and McCann found the same eigenvalues A gk for values of 
m = < 1 close to one. These eigenvalues dissolve into continuous spectrum if m 

becomes small. At the threshold m — 1 between fast diffusion and porous medium 
dynamics, after rescaling with a + 1, we recover the spectrum of the Ornstein- 
Uhlenbeck operator —A + z ■ V. 


Proof of Proposition 6.1 The discreteness of the spectrum follows from the fact that 
C has a compact resolvent, see e.g., [32, Ch. Ill, §6.8]. The latter is a consequence of 
standard Hilbert space methods, cf. [47, Appendix]: For every / G L there exists a 
unique w satisfying Cw+w = f. Via the energy estimate ||Vu ;|| L 2 


< 


Ll 


and the compact embedding H^ a+1 C L 2 a , we deduce that (£ + 1) 1 is well-defined 
and compact. 
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The eigenfunctions of % computed in [37J are regular and thus S (H) C E(£). On the 
other hand, by the energy identity for the eigenvalue equation || Vw\\ L 2 + ^ = A||u>||i ,2 , 
every nonconstant eigenfunction of C is also an eigenfunction of Ft. The statement 
of the proposition is thus proved. ■ 


We can deduce the long-time limit of solutions to the perturbation equation (5.1) 


from the known limiting behavior of solutions in the original variables. Since changes 
of mass are allowed after the change of variables, we have to allow for perturbations 
of arbitrary mass. Solutions to the self-similar pressure equation v(t,x) converge to 
a Barenblatt profile of some radius r which is determined by the initial condition, 

v(t,x) ~ -(r — \x\ 2 ) + as t 1, 


cf. (1.5). By the change of coordinates introduced in Section [2] above and made 
rigorous in Section 7] below, it is v(t,x) = ^(1 — |x |) 2 + w(t,z ) + \w(t,z) 2 in the 
positivity set of v. Hence, considering the long-time limit and using that ||w(t)||L°° < 
h < 1 , we deduce that 

w(t, z ) —>• y/r — 1 as t 1 . 

In the following, we will study the higher order asymptotics of this limit via invariant 
manifolds. As outlined in the introduction to the previous section, we will first 


construct invariant manifolds for the truncated equation (5.2) and then carry the 


results over to the original equation with the help of smoothing estimates. 

The formulation of the invariant manifold theorem requires some preparations. The¬ 


orem 5.1 guarantees the existence of a semi-flow 


Qt 


Li^Li 


defined by S t e5 (g ) = w(t), if w denotes the solution to the truncated equation (5.2) 
with initial datum g. The map [0, oo) x L 2 a 3 ( t,g ) <->• Sl s (g) G L 2 is continuous. 
Following Koch’s (among others) approach to invariant manifolds based on discrete 
semi-flows, see [35], it is enough to study the corresponding time-one maps. That 
is, we consider the mapping S = Si s : L 2 —> L 2 . By the construction in Theorem 


5.1, it is clear that S is Lipschitz, 


ll*%i) - S(g 2 )\\ L 2 < ||gi - g 2 \\i£, 

for any two gi, g 2 G L 2 . Moreover, the following smoothing property holds: 


( 6 . 2 ) 


Lemma 6.1. There exists e 0 > 0 such that for every 0 < e < £ 0 an d 0 < 5 < 1 with 
y/2 (e + h) < 1 the following holds: For any g 1 , g 2 G L 2 fl B £0t s, it holds that 

ll>S'(5 , i) - ‘S'(5 r 2)||co,i < ||gi - < 72 1|■ 


Proof. Let us start by choosing e 0 , e and S as in Lemma 5H. We will prove the 
stronger statement 


\d*dd(w 1 (t,z)-w 2 (t,z))\ < Wgi-teWi*, 


(6.3) 
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for all (t,z) G (|, l) x -£>i(0) and k G N 0 , /3 G N^, where tci and tt7 2 denote the 
solutions to the nonlinear equation (5.1) with initial values g\ and < 72 , respectively. 


We argue similarly as in the proofs of Lemmas 4.8 and 4.9 
from the estimate 


That is, we deduce (6.3) 


l|d t *9f(<Oi - < llsi - Salltj. 


(6.4) 


For this, consider a smooth temporal cut-off function 77 such that 77 = 1 in [|, l] 
and 77 = 0 in [0, . Then 77(7771 — w 2 ) satishes the equation 

dt (77(7771 - w 2 )) + c (77(7771 - 777 2 )) = 77(/1 - f 2 ) + d t rj(wi - tt7 2 ), 


where _/) = /(i< 7 j, Vrc*) for i = 1, 2. By Lemma 4.6, it holds that 

IIV(77 (777i - 777 2 ))|| jL 2 (L 2) + ||pV 2 (77(777i - 777 2 ))|| i 2 (L 2 ) 
< Mfl - /2)||l2(L2) + ||777i - 777 2 || L 2 (L 2). 


Thanks to the decay estimate (2.6) in Theorem 2.1, it holds that 

\\d%d e q d;F(w, V777)|| Loo ^i ^ xBi(0) ) < £, 

for any t G N 0 , (3, 7 G and any solution w G B £ .s- We thus have 


(6.5) 


lh(/l - /2 )||l 2 (L2) 

< ^|| v(77(7771 - 777 2 ))|| L 2 (i 2) + £||pV 2 (77(7771 - 777 2 ))|| L 2 (L 2) + ||777i - 77J 2 || L 2 (i 2). 

Because ||tt7i—tt 7 2 ||z,°°(z.2) < llffi - S^II l 2 holds as a by-product of the proof of Theorem 


5.1. the latter implies (6.4) for (k, |/3|) = (0,1), provided that £ is sufficiently small. 


Moreover, since 


(v(fl - /s)) IU 2 (L 2 ) < £^2 \\d2 - 777 2 )) ||l 2 (L 2 )) 

1<0 


we can argue as in Lemma 4.8 (distinguishing “angular” and “radial” derivatives) 
and derive (6.4) for k = 0. The control over the temporal derivatives is obtained 
by using the equation. Finally, (6.3) follows from (6.4) by the use of Morrey-type 
estimates. This concludes the proof. ■ 

The following decomposition simplifies the discussion. If we denote by L — e~ c : 
L 2 —y L 2 the time-one map corresponding to the heat semi-group, then we write 

S = L + R. 


It is clear that R : L 2 —y L 2 is Lipschitz because both S and L are. Furthermore, R 
is a contraction on L 2 provided that e is chosen sufficiently small and differentiable 
as a map from C * 0,1 to L 2 and from C * 0,1 to C 0,1 : 
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Lemma 6.2. There exists £q > 0 such that for all 0 < £ < £q and 5 > 0 with 
y/2 {e + 5) < 1 the following holds: For any gi, g 2 G L 2 a , 


\\R(gi) - R{ 92 )\\ia < e\\gi - 92 Wli- 


Moreover, 


for all g G C 0 ' 1 . 


imhi <\\9\\l ip and \\R(g)\\coa < hfcon 


The contraction property will be a crucial ingredient in the construction of the 
invariant manifolds below. We will use the quadratic bound on R to show that the 
center manifold meets the eigenspaces tangentially at the origin. 


Proof. In the beginning, we choose £q, £ and 5 as in Lemma 5.1 We first prove the 
Lipschitz estimate. Let w denote the solution to the equation 

d t u> - p _cr V ■ ( p a+1 \7iv ) = f £>s (w 2 ) - f e , s (u>i) 

with zero initial datum. Here ■uy(f) = Sl s (gi) and therefore tc(l) = R(g2) — R(gi)- 
Using the elementary estimate \F £iS (qi,pi) - F £)5 (g 2 ,P 2 )| < £ (\qi ~ + |piP 2 I), 

we compute 

Jw(l) 2 dp a < J'(F etS (wi,Vwi) - F £:S (w 2 ,Vw 2 )) 2 dg u+1 dt 

< £ 2 J J(w 2 — Wi) 2 dg a dt + £ 2 J j \Vw 2 — Vwi| 2 dg a+1 dt. 

The fixed point argument, in the construction of solutions to the truncated equation, 


5.1 


Theorem 


first statement. 


yields that the latter is bounded by £ 2 \\gi — g 2 \\ 2 l2 - This proves the 


If, however, gi = g and g 2 = 0, then the previous argument also shows that 


w(l ) 2 dg a 


< 



F(w, Via ) 2 dpcr+idt, 


with w(t) = St s (g). Since \F St s(q,p)\ < \p\ 2 , we deduce the first quadratic bound on 
R via Lemma 15.11 


Moreover, because fu(0) = 0, Theorem 4.1 yields the estimate 

IMkip ^ WfsA w )\\Y(p), 


where p as in the assumption of the theorem. By Lemma 5.1, the solution w(t) = 
Sl s (g ) remains unchanged if we lower e (if necessary) such that £ < 11 p 11 Lip - Now, 
applying estimates (5.5) and (5.6[) from the proof of Lemma 5.1, we further have 


\\fs,s(w)\\Y( P ) 


< 


Lip 


< 


11 fl 1 11 Lip) 
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provided that e is sufficiently small. We thus have 

ll^)llLip<NlL P - 

Thanks to the elementary estimate ||CIU°° < IICIIid + HCIkip, the second quadratic 
bound follows. 


We are finally in the position to begin our dynamical systems argument. 

The operators £ and L share the same eigenfunctions and the spectrum of L is 
given by the eigenvalues e~ Xek , with the A^’s given in Proposition 6T above. We 
relabel these eigenvalues by {Aj}., g N 0 , where X 3 < A J+ i, that is, eigenvalues will not 
be repeated. The eigenfunctions of £ (and thus of L) can be chosen in such a way 
that they form an orthonormal basis of L 2 . By construction, they are orthogonal in 

Hl + 1 - 


From now on, we keep K £ No arbitrarily fixed. We denote by P c the spectral pro¬ 
jection onto the (finite dimensional) subspace of L 2 spanned by the eigenfunctions 
corresponding to the eigenvalues {Ao,..., Ax}- Furthermore, we set P s I — P c 
and define E c = P C L 2 and E s = P S L 2 , so that L 2 = E C (B E s . We decompose L ac¬ 
cordingly by setting L c = P C LP C and L s = P S LP S . The subscripts “c” and “s” stand 
for “ center 3 ' and “ stable ”, respectively. It is clear that L c has a bounded inverse and 


|| L -€|| < e t\ K for all £ e Nj 


where || • || denotes the operator norm. Indeed, if 11^. projects onto the eigenspace 
spanned by the kt\\ eigenfunction then U^wit) = e~ Xkt Ilkg- Hence, ||L“ M\h = 

\\Pc9W% = Efc<x( n W,'0fc)l2 = Efc<jf e2Afc ( n fc u; >V’fc)ia < e 2Ax ||u ’||^2 ■ The above 
bound follows by iteration. Likewise, 

\\L e s \\ < e~ aK+1 for all £ £ N. 


The latter follows from the energy decay rate for the linear equation and a spec¬ 
tral gap estimate on E s . Indeed, on the one hand d t w + Civ = 0 implies that 
11 P s iv 11 2 jj 2 + ||VF s tc|| 2 2 = 0. On the other hand, Ax+i||w ||^2 < ||Vu)||^ 2 for all 

w £ E s . Combining both inequalities and iteration yields the desired bound. 
Following the notation introduced in [35], we finally define 

A s = e~ XK+1 , A c = e~ XK , and A max = 1. 


With these estimates, we have 

j|£; 1 || < A" 1 , \\L S \\ < A s , ||£|| < A max . (6.6) 

Furthermore, if £ gap G (0,1) is chosen small enough, we can find A_ < A + such that 

A s H“ £gap ^ A_ A c £gap? 

A m ax H“ ^gap ^ A_|_. 
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For any function w E L 2 , we define ||w|| = max {|P c w|U|, ||P a iu||L 2 }. Moreover, for 
any sequence {wk}ke z we set 

||{wfc}fcez||A_,A+ := sup max{A“ fc ||w fc ||, A*||w_ fc ||} , 
fceN 0 

and for any sequence {wfe}fc e N 0 , 

||{u’fc}fceN 0 ||a_,+ : = sup Al fc ||'iUfc||. 

fceN 0 

The corresponding Banach spaces will be denotes by I\_,a + and £a_,+, respectively. 
In a first step, we construct the center manifold. 

Proposition 6.2 (Center Manifold). Let e be small enough so that Lip(P) < £ gap . 
There exists a map 9 : E c —>■ E s with 9(0) = 0 and D9( 0) = 0 such that the 
submanifold 

XV C = {w c + 9(w c ) : w c E E c } 
has the following properties: 

1. (L 2 regularity) The map 9 is Lipschitz with Lip(d) < e gap . Moreover, 

\\0(9c)\\li < hcWli ( 6 . 7 ) 


for all g c E B £ j fl E c . 

2. (C 0,1 regularity) The map 9 from C 0,1 to C 0,1 is Lipschitz on B £ j. Moreover, 


ll^(fl'c) ||c°’i < II 9c 


2 

C 0 ’ 1 


( 6 . 8 ) 


for all g c E B £: § fl E c . 

3. (Invariance) The restriction on W c of the semi-flow {S'* , 5 }t>o can be extended 
to a Lipschitz flow on W c . In particular, Si S (W C ) = XV c for all t > 0, and 
for any g E W c there exists a negative semi-orbit {?/;(£) } t < 0 in W c such that 
w( 0) = g. 

4■ (Lyapunov exponent) A point g belongs to W c if and only if there exists a flow 
{w(t )} te r with u;(0) = g and 


w(t)\\ L 2 


< 

rs_/ 


A+ for all t > 0, 
A(_ for all t < 0. 


5. (Optimal Lyapunov exponent) A point g belongs to XX) if and only if there 
exists a unique flow {w(t)} te n with w(0 ) = g and 


\w(t )\\li 


< 

r^j 


(A max + e gap f for all t > 0, 
(A c - Sgap) 4 for all t < 0. 
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In the statement of the proposition, DO denotes the Frechet derivative of 6. 
The manifold W c is not unique but depends on the cut-off function r) £ 

With slightly more effort, we could prove (6.7) for all g G E c . 


Proof. Most parts of the theorem are an immediate consequence of Theorem 
2.3]. We thus only sketch the proof. The results which are new or adapted to our 
particular study, however, will receive a detailed presentation. 

For any g c G E c and {w e }e £ z we define {J k (g c , {we}iez)}kez by 


Jk ( g c , {u>i}eez) 


P s S(w k - 1 ) + Lf 1 P c (w k+1 - R(w k )) for k < -1, 
P s S(w- 1 ) + g c for k = 0, 

S(w k ~ i) for k > 1. 


It is easily checked that J = {J k } k &z maps E c x £a_,a + to f , A_,A+, and for every 
fixed g c G E c , J(g c , •) is a contraction on ^a_,a+- Thus, by Banach’s fixed point 
theorem, for every g c G E c , J(g c , ■ ) has a unique fixed point {w k } k£ z hi ^a_,a+- It 
is a solution to the discrete semi-flow with P c w 0 = g c and extends to a global flow, 
that is 

w k = S k (w 0 ) for all k G Z and P c w 0 = g c . 

We then set 0(g c ) = {w k } ke z and 0(g c ) = P s 0 o (g c ) = P s w o, that is the projection 
of first coordinate of the sequence {w k } k &z onto E s . Thus 0 : E c —>■ f , A_,A + and 
0 : E c —y E s , and both maps are Lipschitz as a consequence of the fixed point 
argument. Moreover, it can be shown that Lip (0) < e gap . Indeed, using (6.6) and 


Lemma 6.2, via iteration, we observe that 


ll%c)-»&)IU;< | 


+ £ 


k 

g a P XI 
1=1 


A, 




/) | \\o(g c ) - e(~g c )\\A-A + - 


Then the statement follows from using the Lipschitz estimate for 0 and letting k f oo. 
The argument for (6.7) proceeds similar. Notice that \\R(w _ k )\\ L 2 < ||S'i)||Li P ~ 


|tc_fc_i ||22 via Lemmas 


6.2 


and 


6.1 


Hence, via iteration 


li«(9o)iUs< |Lyii%)|| A _, A+ + c]r(LY 1 ||9( 9 c)||L,a, 


e=i 


and the statement follows from letting k f oo via the Lipschitz estimate for 6 


because d(0) = 0. Notice that ( |6.7[ ) entails that D6( 0) = 0, where DO is the Frechet 
derivative with respect to L 2 a . 

That 0 is also a contraction as a mapping from C* 0,1 to C 0,1 can be seen as follows: 
By the definition of 0 and the regularity estimate from Lemma 6T, for any g c and 
~g c in E c , it holds that 

ll%c) - 0(g c )||cm ^ ll^(w-i) - 5(w_i)|| c o.i < ||w_i - u)_i|| L 2 , 
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where {w k }ke z = 9(g c ) and {w k } k &z = 9(g c ). Hence 

II %c) - 9(gf) He 0 ’ 1 < II 0(9c) ~ 0(~g c )h.,A + ^ he - 9c\\u>, 


and the last inequality is due to the Lipschitz dependence of 9 on g c . The L\ norm 
being trivially estimated by the U 0,1 norm, we conclude the Lipschitz regularity with 
respect to C 0,1 . The quadratic bound (6.8) is established as follows: Via triangle 
inequality it holds ||0(<7c)||c'°.i < ||L s u;_i||cm + ||i?(iy_i)||co,i. We use the quadratic 
estimate for R in Lemma (L2 and the smoothing estimates for L and S in Lemmas 
to deduce ||6 , ((yf c )||co.i < ||P s tc_ 2 ||L| + ||xt?— 211 ^, 2 , and the statement follows 


4.9 


and 


6.1 


via iteration as above. 

The remaining parts follow by construction, see 


Proposition 6.3 (Stable Manifold). Let e be small enough so that Lip(i?) < e gap . 
There is a continuous map u : E s xLl —>■ E c such that for each g G W c , u(—P s g,g) = 
—P c g and the manifold 


M g — g + {u(w a , g)+w s : w s e E s } 


has the following properties: 


1. (L 2 a regularity) For every fixed g G L 2 , the map u( ■, g) : E s —>■ E c is Lipschitz 
continuous with Lip(z/( •, g)) < £ gap . 

2. (C 0,1 regularity) The map v is continuous as a mapping from C 0,1 x C 0,1 to 
C 0,1 . Moreover, for every fixed g G C 0,1 , the map u( ■, g) : C' 0,1 D E s —> C 0,1 is 
Lipschitz with Lip■, g)) < e gap . 

3. (Invariant Foliation) It holds Sl s (M g ) C M s t ^ for all t > 0 and 

Mg = l~g g Ll : sup Az k \\S^{g) - < (1 - k) _1 II P s (g -~g)\\ A, 

l fceN 0 J 


where 


k = max 


A s + e gap A_ + e. 


gap 


A_ 


A, 


4- (Completeness) For every g G L 2 a the set M g D W c is exactly a single point. In 
particular, {M g } geL 2 is a foliation of Lf overW c . 


Proof. Up to the proof of the second item and the precise decay rate in the third 
item, the statement can be found, e.g., in [35, Theorem 2.4], For the new results 
(which are peculiar to the situation at hand), we have to review the construction 
of v. For this purpose, consider the mapping / = {I k }ke n 0 : E s x ^a_,+ —; > ^A_,+ 
defined by 


h(g s , {wfjteNo) 


g s +L c 1 P c (wi - R(w 0 )) for k = 0, 

P s S(w k - 1 ) + L~ 1 P c (w k+ 1 - R{w k )) for k > 1. 
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Let g s be fixed for a moment. Then I(g s , ■ ) is a contraction on £a_,+ with contraction 
constant 

A$ T ^gap A_ T ^gap 1 

A- ’ A J ' 

We let R > (1 — k) - 1 ||s , s||l 2 . Since I k (g s + P s g, {S e (g)} eG ^ 0 ) - S k (g ) = sAo for 
every (7 £ L 2 , we notice that I(g s + P s g, •) maps the ^A_ j+ -ball of radius R around 
{S k (g)}kef<io on fo itself. Hence, Banach’s fixed point theorem yields the existence 
of a unique sequence {w k } keNo such that {w k } k£ ^ 0 = I{g s + P s g,{wt } leNo ) and 
satisfying 

||{wfc}fc e N 0 — {S k (g)} ke - No ||a_,+ < R, 

By construction, {w k } ke n 0 is a discrete semi-flow with P s w 0 — g s + P s g. We now 
define u{g s ,g) = {u> fc } fc eN 0 and z/(g s ,g) ■— P c (ta 0 — d)- Via iteration we obtain 


k = max 


v{ds,g) - v(g s ,g )\\ L 2 < 



- Hgs,g) IU_,+, 


and since z>( •, < 7 ) is Lipschitz by construction, we see that Lip(V( •, g)) < e gap by 
letting k t 00 . 

Then, for any two g s and g s in L 2 and {w k } keNo = u(g s ,g) and {u> k }ke n 0 = v(9s,g), 
it holds that 


I "(gs, g) - v(g a , g)\\c°'i = \\P c (w 0 - w 0 )\\ c oa < \\L c P c (w 0 - w 0 )\\ c ^, 


by the invertibility of L c in C 0,1 . We now apply Lemma 4.9 to deduce 


I W(g s ,g) -v{g s ,g)\\c°’i % H^cOo - w 0 )\\ L 2 = \\u(g s ,g) -u(g a ,g) || L a. 

Invoking the Lipschitz bound on ■ , g) in L 2 and the trivial embedding of C 0,1 into 
L 2 , - ,g) is Lipschitz with constant < £ gap . A very similar argument shows that 

v is also C 0,1 Lipschitz in the second argument. 

Based on the construction of {w k }ke n 0 , most of the properties in the statement 
of the proposition are readily verified. For instance, invariance follows from the 
construction by uniqueness. Moreover, if g — g + u(g s , g) + g s E M g for some g s £ E s , 
then g = wq where {rUfc}fc 6 N 0 = v(g s , g). In particular w k = S k (g ) for all k £ N 0 and 
thus ||{S' fc (fi , )}fceN 0 - {S k (g)}k£N 0 ||a_,+ < R, with R = (1 - k) _1 | \P a (g - g)\\ l%- On 
the other hand, if the latter holds for some g £ L 2 , then a similar argument shows 
that g £ M g . 

Finally, completeness can be seen as follows: Let x '■ E s x L 2 —» E s be given by 


X(g s ,g) = 0(u(g s - P s g,g ) + P c g). 

Since 9 and u( ■ , g) are both contractions with respect to the L 2 topology, so is 
x( ■, g)- Hence, for every g £ L 2 , there exists a unique fixed point g s = x(9s, g)- We 
now set 

9c ■= v(g s - P s g , 9 ) + P c g, and g s = 9(g c ). 
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In particular, the definitions imply that g c + g s £ W c fl M g . However, every intersec¬ 
tion point is a fixed point of x (' , <?)■ Hence, by the uniqueness of the fixed point, 
W c fl M g = {g c + 


In the previous propositions, we have derived an invariant manifold theorem for 


the flow of the truncated equation (5.2). This theorem can be carried over to the 
original problem. 


Proof of Theorem\2.2. That W x ° c is locally flow invariant is clear from Theorem 


2.1 


and Proposition 6.2 For the same reason, Wl oc contains all the negative semi-orbits 


that satisfy the constraint. 

Let g be the unique point in M g fl W c found in Proposition 6.3 We claim that 

SU'g) = S‘(g) for all t > 0 , 


(6.9) 


if 0 < £ < £* for some £* sufficiently small if | < 7 11 Lip < e 0 for some 0 < £ 0 < e 
sufficiently small. From this, the characterization of the stable manifold M g in 


Proposition |6.3| yields that 
11^(0) -S\g)\\Ll < max 




Ar 


A_ - A, - £ gap ’ A r - A_ - e. 


gap 


A f _ for all t > 0. 


We now invoke the smoothing estimate from Lemma 6 H to deduce the statement of 
the theorem. 


It remains to verify (6.9). The construction of W c and M g in Propositions 6.2 and 


6.3| yields the existence of points g c G E c and g s G E s such that g = g c + g s and 
g c = v(g s — Psg, g) + P c g and g s = 9(g c ). Hence g c solves the fixed point equation 


9c = v{0(gc) ~ P*g,g) + P c g = fy% c ) - P s g,g) - v{o,g). 


( 6 . 10 ) 


With respect to the C 0,1 topology, 6 is Lipschitz in B s $ by Proposition 6.2 and 


v(-,g) is a contraction by Proposition 6.3, provided that e is chosen sufficiently 


small. It follows that the mapping g c —y o(6(g c ) — P s g , g) + P c g is a contraction with 


respect to C 0,1 in B St s, and thus, in B Si s, ( 6 . 10 ) has a unique solution g c = g c (g) 
which depends continuously on g by the C 0,1 continuity of u. Now, since 0(0) = 0, 
u( 0, 0) = 0 and ,r/ c (0) = 0, by continuity, there exist 5 > 0 and s > 0 such that 


g G B s j if g G B-g. This proves ( |6.9[ ). 
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7 Applications of the invariant manifold theorem 


In this final section, we give four applications of the invariant manifold theorem. 
We prove the long-time asymptotics for solutions to the perturbation equation with 


small initial data: solutions to (5.1) converge to (small) constants, see Theorem 2.3 
above. We refine the estimate that is deduced from Theorem 2.2 and compute the 
precise rate of convergence. On the level of the porous medium equation, this rate 
is saturated by spatial translations of the Barenblatt solution. This first result is 
obtained by choosing K = 0 in the notation of the previous section, that is, we 
investigate the asymptotics towards the one-dimensional manifold spanned by the 
constant eigenfunctions corresponding to the eigenvalue Ao = 0. In the second exam¬ 
ple, we fix the long-time limit of solutions by letting the mean vanish asymptotically. 
Hence, the long-time limit is the zero function. According to the invariant manifold 
theorem, solutions converge to the 1 + N dimensional center manifold corresponding 
to the eigenvalues A 0 and \\ (thus K = 1) with an improved rate A 2 — £, see The¬ 
orem IH On the level of the porous medium equation, we subtract out the spatial 
translations. Climbing the eigenvalue ladder one rung up (K = 2), in Theorem 
we additionally asymptotically suppress the center of mass of solutions, so that 


2.5 


the second order corrections (affine transformations) in the asymptotic expansion 


becomes accessible. Finally, in Theorem 2.6, we consider the dynamics on the rate 
A 4 (hence K = 3) in the case where N < <7 + 1 , which, on the porous medium level, 
is governed by dilations. 

In principle, we can access all eigenmodes through the invariant manifolds of Theo- 
Since the mathematical arguments behind the computation of even higher 


rem 


2.2 


modes do not change significantly, we will finish this study with the third order 
corrections. 


The proofs of Theorems 2.3 to 2.6 are in order. 


Proof of Theorem 2.S\ With regard to the notation introduced in the previous sec¬ 
tion, we suppose that K = 0, and thus E c = span{'0 olo } = R. We choose A between 
Ai and A 0 and find £ gap > 0 such that A_ = e~ x e (e~( g + L + e gap , 1 — £ gap ). Let 
e, 6, £ and 6 be given as in the statement of Theorem 


2.2 


Our first claim is 


IF r loc 2* {a E R : |a| < 5} . 


(7.1) 


Indeed, on the one hand, if w G H / c loc , then w = a for some |a| < 5 , which shows 
the inclusion “c”. On the other hand, if |a| < 6 and w(t) = 6 for all t < 0, then 
{tc(t)}t<o is a negative semi-orbit meeting the norm constraints. Hence, by Theorem 
2.2, w(t) E for all t < 0. This implies the inclusion “d”. 

In particular, by Theorem 2.2, there exists a unique semi-flow {to(t)} t >o in W l ° c 
such that 


| w(t) — w(t) He 0 ,! < Al for all t > 0. 


(7.2) 


However, (7.1) implies that w is constant with modulus less than S, say w = a. 
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It remains to optimize the rate of convergence to e ( <T+1 )*. Consider now a(t) = 
I( 0 ) 1^ 1 1 w (t) A short computation reveals that 

d f 

—a(t) = P / F(w,'Vw)dp a+1 , 


and thus, taking into account F(q,p) < |p| 2 and the above decay estimate (7.2), we 
deduce that 0 < < e~ 2Xt . Hence, there exists a constant a G R such that 

\a(t) — a| < e _2At . On the other hand, | a(t) — a\ < e~ xt by the definition of a(t) and 
(7.2), and thus d = a. 

The above computation shows that P s w(t ) = w(t) — a(t) G E s satisfies 

WPgwit) ||c°’i e ~ Xt f° r all t > 0. (7.3) 

(Notice that E s is the space of all mean-zero L 2 functions.) We shall go over to the 
discrete semi-flow for a moment, that is, we set Wk = w(k) for k G No- Using the 
decomposition S = L + R introduced in the previous section, we write 


w k = L k w 0 + J2 Lk ~ jR ( w j-i)- 

3 = 1 


6.2 


From a variant of Lemma 
which applied to the above identity, yields 


we deduce the estimate \\R(g )\\ L 2 < £ 1 P\\P s g\\l[l\\P s g\\ L 2 , 


\P s w k \\ L 2 < A k s \\P s w 0 \\ Ll +C'£ 1 / 2 ^A^||P s u; i _ 1 ||i/ p 2 ||P^_i|U 2 

3 = 1 


(7.4) 


for some C > 0. We claim that for e small (but independent of k) it holds that 

II-Ps'UT-IIl 2 < 2CAj||t(; 0 || L 2 for all k G N 0 . (7.5) 


We argue by induction. The statement if trivial for k — 0,1. Assuming that (7.5) 
holds true for i = 1,..., k — 1, we deduce from (7.4) 


\ P sWk\\L% < ( 1 + 2C' 3 / 2 £ 1 / 2 A s 1 \\PsWj-iWul ) A^||wo||l 2 - 

3 = 1 


Because ||P s u;j_i||Lip is exponentially decaying by (7.3), the statement (7.4) is true 
for any k provided that £ is sufficiently small. From (7.5) we infer the statement 
for the continuous flow with the help of smoothing estimates, cf. (6.3) in the proof 
Lemma 16.11 

(The derivation of (7.5) was inspired by a similar argument in the proof of 
Theorem 1.1].) 
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Proof of Theorem 2.f. Suppose that K = 1 and thus E c = span {1, Zi,..., zn}- 
Moreover, A s = e _2 ( CT+1 ) and A c = e~^ a+1 \ We choose A G (cr + 1, 2(a + 1)) and set 
A_ = e~ x . There exists a small e gap > 0 such that A 6 (A s + e gap , A c — £ gap ). Let 
0 < £ < e, 0 < A < A, and g be given as in Theorem 2.2 Then, setting w(t ) = S^fg), 
it holds 


||u;(t) — w(f)He 0 ’ 1 ^ e xt for all f > 0 . 
Because iu(t) G W c , we find a{t ) e R and 6 (f) G R A such that 

w[t) = a[t ) + 6 (f) • z + 8(a(t),b(t), z). 

The statement now follows from the estimates 

\a(t)\ < e -2( ' <T+1 - )t for all f > 0 , 
| 6 (f) -e~ ia+1)t b\ < e - 2(ff+1)t for all f > 0 , 


(7.6) 


(7.7) 

(7.8) 


and the quadratic bound for 6 in (6.8h 


The proof of (7.7) proceeds in the same way as we obtained the estimate for a(f) in 


the proof of Theorem 2W Notice that 

1 


a(f) = 


\Bi 


w(t , z ) dn a (z) —y 0 as f f oo 


by the virtue of Theorem 2.3 and the assumed long-time limit of w(t). Now, since 


w solves the perturbation equation, it holds that 

[ F(w,A/w)(t)dfi a+ i < I |Vw)(t)| 2 d/v+i, 


d ~(+\ P 
Jt a(t) = w 


and thus, using the decay rate for \7w in Theorem 2.3, we deduce (7.7) 


The second estimate (7.8) is established in a similar way. We notice first that 


6 j(f) = c a , N J wzi dg a , 

for all i G {1,..., N} and some constant c a and thus, differentiation and integra¬ 
tion by parts yields 

d o 

= -(cr + l)6j(f) + c a)N (l3 + 1 F(w , Vw)Zi dHa+ 1 - 


With the help of Theorem 2.3, the latter leads to the estimate 

d 


j t ( e( " +1) ‘bo) 


< p-O+W 


and we deduce (7.8) with b = lim^oo e^ +r>t b{t). This concludes the proof. 
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Proof of Theorem \2.5[ As for the notation of the previous section, we are now in 
the situation K = 2. Compared to the previous theorem, we augment the finite 
dimensional projection space E c by the eigenspace associated to the new eigenvalue 
A 20 = 2(cr + l). This eigenvalue corresponds to affine transformations of the attract¬ 
ing ground state and the eigenfunctions are harmonic homogeneous polynomials of 
degree two, which can be expressed in the form if(z) = z ■ Az for some symmetric, 
trace-free matrix A. We will further split A = B + C, where B is off-diagoual and 
symmetric and C is diagonal and trace-free. Therefore, if w is as in Theorem |2.2[ it 
holds 

w(t, z) = aft) + bft ) • z + z ■ B(t)z + z ■ Cft)z + 9 (aft ), bft ), Bft), Cft), z), 

for some aft) E R, bft) E R N , Bft ) E R NxN off-diagonal, symmetric and C(t) E 
R NxN ( p a g ona i ; trace-free. Besides showing convergence for Bft) and C(t), we n eed 
new estimates for aft) and bft). Repeating the estimates form Theorem 


2.4 


and 

using the assumptions, we first notice that \a(t)\, | 6 (t)| < e - 2 (o-+i )t s j nce 5 — q c f 


(7.7), (7.8), and thus ||iy(t)|| c o,i < e d* from the statement of Theorem 


2.4 


Let {E n } and {F n } be orthogonal bases for the spaces of off-diagonal symmetric 
matrices and trace-free diagonal matrices, respectively. Then B(t) = 'Yh n Pn(t)E n 
and C(t) = f° r some Pn(t), 7 n (t) E R. It is clear that 


finif) C ( 


N r 

E E z / 

i,j =1 J 


ZiZjwft) dn a , 7 n (t) = c aN 


N 

* 7=1 


ZiZjw(t ) (i/Xo-, 


where c^at may be different in both formulas. Using the perturbation equation and 
integrating by parts a couple of times, we compute 


f ZiZjwdno 


= 2 8 * 


w dfi a+ 1 - 2(a + 1) / ZiZjW dfi a + (/3 + 2) / ZiZjF dfjb a+1 , (7.9) 


where F = F(w, Vtc) < | \7w\ 2 is bounded by e _ ( 4 ( fT+ 1 ) _ 2 d t . First, if i 7 ^ j, then the 
above implies that |B — e 2 ( <T + 1 ) t 5(i)| < e - 2 (o-+i- 6 )t f or some off-diagonal symmetric 
matrix B E R NxN . Similarly, because ^ ■ Ff-dij = tr (F n ) = 0, it holds |C — 
e 2 ( <T+ i) t C'(f)| < e - 2 ( CT + 1 - £ ) t f or some diagonal trace-free matrix C E R NxN . The 
above bounds on B(t) and C{t) reveal that Aft) = B(t) + C(t) is bounded by 
e — 2 (<r+i)i. u s j n g this new estimate in the expansion of iu and applying the quadratic 

2(cr+l )t 


bound on 9 in ( 6 . 8 ) yields that ||w(t)||c?o,i 


< 


As a consequence, we can 


choose e = 0 in the previous computations. We have thus shown that 


A(t) - e~ 2{a+l)t A 


< 4(cr+l)i 


where A = B + C and |a(£)|, \b(t)\ < e - 4 (o-+i)t_ j n particular, the both, the affine 
term and 9 in the expansion of w are of higher order. This proves the statement of 
the theorem. ■ 
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Proof of Theorem \2.6[ We consider the case K = 3, that is, compared to Theorem 
m the new eigenvalue Aoi has multiplicity one and the corresponding eigenfunction 
is V’on (z) = 1 — 'f\z\ 2 , where 7 = IV -1 (2 (cr + 1) + N). Consequently, 

E c = span {l zn, 1 — 7 |z| 2 } U span {7 ■ Az : A = A r , tr(H) = 0}, 


and thus, according to Theorem 2.2, there exists a solution w to the perturbation 
equation which is of the form 

w(t, z ) = a[t) + b(t) ■ z + z ■ A(t)z + c(t) (l — q|z| 2 ) + 6 (d(t), b{t ), Aft), eft), z'j , 


and which satisfies 


\w 


( t ) — u>(f)|| c o,i ^ e At for all t > 0 . 


(7.10) 


By the orthogonality of the eigenbasis, eft) has the representation 

c(t) = c Nt<T J wft) d/ip cjv , o -7 J w(t)\z\ 2 d/j a , 

for some constant c Na > 0. It is clear that |a(t)|, \bft)\, \A(t)\ < e — 4 ( cr + I )* by the 
proof of the previous theorem. Hence, recalling the formula for aft ), it holds that 
eft) = — ctv,o-7 / w(f)|z | 2 d/i 0 - + 0(e _4( ^ cr+1 i f ). To determine the decay behavior of eft ), 
we use (7.9[) in the form of 


^ [ \z\ 2 wdn a — — (2(u + 1) + IV) + f w d/j^ + (f3 + 2) f \z\ 2 F d/j, a , 


where F = Ffw,'Vw). Since H = 0 in the statement of Theorem 2.5, it holds that 
ll^llc 0 ’ 1 e -(s(o-+i)-e)t_ j n particular, since F is essentially quadratic in |Vw|, the 
third term on the right-hand side of the above identity is of higher order. Recalling 
the decay behavior of aft), we thus have 0 < ^ (e( 2 ( tr+1 ) +JV ') t c(i)) < e N ~ 2 F+ l )t m 
Because N < 2(cr + 1), there exists thus a constant c G R such that 

I eft) - e-( 2 (fT+ 1 )+Ar )*c| < e“ 4(CT+1) h 


Hence, by the virtue of ( 6 . 8 ) and the estimates on aft), bft), Aft), the statement of 
the theorem follows from (7.10) and the eigenvalue expansion of w. ■ 
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